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This document is the third in a five volume report which describes a 
comprehensive digital computer simulation of the dynamics of heavy lift 
airships and genetically similar vehicles* 

The work was performed by Systems Technology, Inc*, Hawthorne, 
Call^'ornia for the Aeronautical Systems Branch in the Helicopter and 
Powered Lift Division of the National Aeronautics and Space Administra- 
tion, Ames Research Center, Moffett Field, California. The simulation 
development was carried on between September 1979 and January 1982 and 
is currently instc? 1 led on the Ames Research Center CDC 7600 computer* 
The contract technical monitors for NASA were Dr. Mark Ardema, Mr. Alan 
Faye, and Mr. Peter Talbot. STi's Program Manager was Mr. Irving 
Ashkenas. 

The authors wish to acknowledge the technical contributions of 
Mr. Robert Heffley, Mr. Thomas Myers, and Mr. Samuel Craig and the fur- 
ther contributions of Mr. Allyn Hall, Ms. Natalie Hokama, and Ms. Leslie 
Hokama in simulation software development. Special thanks are due to 
Ms. Kay Wade, Ms. Linda Huffman, Mr. Charles Reaber, and STI's produc- 
tion depactn^nt for the preparation of the five volumes of this report. 
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The heavy lift airship alaulation consists of three coaputer pro- 
grams* The vehicle Is aodeled as a buoyant envelope (hull) with lift 
propulsion units (LPUs) attachcJ. The LPUs each consist of one rotor 
for upward lift* one propeller for forward propulaloni and a fuselage 
(nacelle). The aodel also Includes landing gears, an attachable slung- 
payload, and a mooring point* The characteristics of ail aodel eleaents 
are user defined* 

The three programs are: 

HLASIM — Models the powered vehicle In flight* 

HLAMOR — Models the moored unpowered vehicle con- 
strained to a mast at the mooring point* 

HLAPAY — This program la the largest of the three* 

It models the vehicle In flight carrying a 
slung payload* 

This simulation was divided into separate programs because: 

1) This allows modular development and testing of 
the program during development and future altera- 
tions* 

2) T,^,rge parts of the main program were unnecessary 
for program HLAMOR because it Is unpowered, and 
has only three degrees of freedom* The trim for 
HLAMOR is consequently completely different* 

3) The payload was developed and tested as a separ- 
ate program. HLAPAY executes the payload program 
in addition to program, HLASIM* Having a vehicle 
only program (HLASIM) frees the user from having 
to prepare payload Input files unless the vehicle 
with payload is to be modeled* 

All programs have the same basic algorithm and large parts of the 
programs use the same code* They also share the same datr files when- 
ever possible. 
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The basic program algorithm is: 

1) Read and initialize data* 

2) Calculate the trim conditions. 

3) Output the results of the trim calculations* 

4} Calculate and output the stability derivatives if 
requested* 

5) Repeat 1'4 for additional trim conditions If 
requested* 

6) Read and Initialize the time history data. 

7) Write the output data at time 0.0* 

8) Integrate the vehicle states forward in timet 
implementing appropriate control commands and 
gust disturbances* 

9) Write program output data at the user requested 
print intervals. 

10) Repeat 8 and 9 until the end of the time history* 

A. KSCRIPTKNI OF THIS AMD CTHBR BELATED MRNOALS 

An overview of the various manuals associated with this simulation 
is in order before continuing with detailed discussion of the program. 
In addition to this User's Manual, there is a Programmer's Manual and a 
Technical Manual* 

This User's Manual is designed to provide the user with the basic 
information necessary to run the program as it has been designed. This 
manual does not vUscuss any of the internal workings of the code or the 
technical details of the equations and their derivations* This manual 
describes the various data files necessary for the program and explains 
the output from the program and the various options available to the 
user when executing the program. The discussion of the data files is 
limited to: 

1) The type of data contained in each file* 

2) The Inputs necessary to create special conflgura** 
tions. 
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3) Inputs vhose nature la spaclalliad or not obvl* 
oua« 

4) Additional data file inforaation is contained in: 

Appendix A - tabulates all input variabler;* 

It indicates which values are valid for the 
various variables and other special considera- 
tions which the varlablee nsy have* 

Appendices B and C - contain sanple sets of 
input files and the output resulting froa 
those input files* 

l>e Prograaner*8 Manual is designed for the ■aintenance prograMer 
who will uc supporting this program. It explains the logic of the vari- 
ous program im>dules, and in some cases gives a detailed explanation of 
the reasons ft.: various Implementations. The topics discussed in the 
User *8 Manual will not be repeated in ihe Programmer's Manual * Conse- 
quently^ the maintenance programmer will have to consult both of these 
manuals when working with the program* The Programmer's HaQual conteis? 
several appendices^ including a dictionary of program varlables» a list 
of all subroutines and their purposes » a subroutlne/coismon block/crosa 
reference listing, and a calllng/called subroutine cross reference 
listing* 

The Technical Manual contains a detailed discussion of all simula- 
tion models, Including derivations of all the equations, and methodology 
for calculating the required program input data. The user will have to 
consult tikis manual for all technical information he requires in gener- 
ating the input data files or in understanding the output. 

B. NOTES TO A FIRST TIHE OSER 

The heavy-lift airship simulation programs are rather complex, 
containing many user controlled options. Therefore, it is recommended 
that any configuration to be simulated should first be analyzed using 
program HLASIM. By using this program, the user can def.'ne geometry and 
mass characteristics and make test runs to validate his input files* 
The biggest problem confronting a new user will be the generation of 
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data files which accurately represent the vehicle he wishes to aodel* 
Several runs should be oade to check the output for reasonable results. 

Maintaining data files can be a significant problen. To assist in 
overcoBing this bookkeeping task, input routines were designed to allow 
the sane data files to be used by the three different prograas. For 
example, the mooring program (HLA)K}R), which does not use the flight 
control system inputs, will still read the same data file (HISDTA) as 
program HLASIM. HLAMOR passes over the unwanted data. In this way, the 
user can analyze the flight conditions and the moored conditions of the 
same vehicle, confident that the configuration input data was exactly 
the same for both prograas. 

It is recommended that the user create a basic set of input files. 
Then, make copies of those files and alter the data appropriate to the 
program run wanted. These alterations can be noted in the comment lines 
of the interactive questions data file, and later the run can be easily 
identified by knowing chat the basic data with Che noted changes was 
used. This method of copying files Insures that all runs start with the 
same basic data. If the user alters the basic data file he runs the 
risk of falling to return the data to its original form and thereby 
invalidating future runs. 

The following section describes the minimum data files necessary to 
execute program HLASIM. This will allow a new user to run the program 
with minimal initial effort. Additional data files required for special 
options, can be addod as they are needed. 

Th»> following eight data files are required to execute the program 
with a minimum of options.* 

The first two files define the desired printed output varlsbles and 
the basic program user options: 


*In order to execute the program with this minimum data set, the 
user must enter "False" to the Interactive question "Time History?" 
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1) OUTLST 

2) Interactive responses (INPUT) 

The next file which is supplied along with the program code, con~ 
tains the program error messagetc 

3) ERMSSG 

The remaining five files are used to define the basic vehicle con- 
figuration and trim conditions: 

4) GMOTA 

5) ARODTA 

6) INFDTA 

7) PLMDTA 

8) TRMDTA 

Using these files, the program will calculate the trim and print the 
results. 

Additional Data: 

9) HISDTA — Time history data file. This file must 
be supplied if the user has answered the inter- 
active question "Time History?" as true. 

Gust string files (RGl, RG2, RG3, and RG4) are not necessary as long 
as the user Inputs GSTFLG ■ F in the data file HISDTA. 

This manual describes each of the above files in Section 2. The 
user should consult these sections along with Appendix A, B, and C when 
generating data files. Section 11 of this manual contains the specifics 
of how to create and load the Job to run the program on the NASA/Ames 
CDC 7600/Scope System. 
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section discusses each Input data file necessary to the program 
(HLASIM)* The discussion Is not exhaustive* The Methods for estlMtlng 
<:hc engineering constants are found In the "Technical Manual" and the 
definitions of each Input variable Is given In Appendix A of this nen*- 
ual* Sample data files are given In ^pendlx B* The discussion here 
Include the file foraat and special or not obvious uses nade of any var- 
iables. 

The follovlng Is a complete list of the Input data files: 


GMDTA 

ARODTA 

IFCDTA 

PLMDTA 

TRMDTA 

HISOTA 


Vehicle* Trim* and Time history 
variables (Discussed In Section 2-A 
to 2-F) 


MORDTA 


Mooring exclusive data (discussed In 
Section 8-B 


PAYDTA 


Payload exclusive data (discussed In 
Section 7-B) 


RG1-RG4 Gust strings 

RGS, RG6 Guot strings (discussed In Section 2-G) 


OUTLST 


User selected output variables code 
numbers (Vehicle discussed In Section 2-H) 


PYOUTL 


User selected output variables code 
numbers (Payload-discussed In Section 7-B) 


INPUT 


Run Parameters (discussed In Section 2-1) 


ERMSSG 


Error messages (Supplied with program discussed 
In Section 2-J) 



( 
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A. DATA PILE: GMDTA 


Contents: Geometry and Mass Characteristics 

The variables In this file are listed under the following output 
headings : 

e "GEOMETRY INPUTS," 
e "MOORING POINT GEOMETRY," 

e "LANDING GEAR ATTACH POINTS AND SPRING CONSTANTS," 

e "MASS AND MOMENT OF INERTIA INPUTS" 

The definitions and descriptions for these Input variables are found In 
Appendix A. The moment of Inertia characteristics are discussed In the 
technical manual accompanying this program. Consequently, only those 
Inputs which have special significance or are used in a special manner 
are discussed in this section. All Inputs are in Namelist format. 

Figures 1 and 2 Illustrate the geometry configuration vectors and 


their reference centers s 

The legend 

for Figs. 1 and 

2 is as follows: 

?i - RTALOC 

GMDTA 

ry - RACELC 

HISDTA 

?2 - RATCH2 

GMDTA 

rg - RVSNLC 

HISDTA 

r3 - RLTCH2 

»1DTA 

?9 - RR0TR2 

GMDTA 

?4 - RMORPT 

GMDTA 

rio - RPR0P2 

G>1DTA 

rs - RATHGl 

GMDTA 

rn - GGLP2 

GMDTA 

rg - RHULCG 

GMDTA 

ri2 - RACLP2 

AROOTA 


P^ - Reference Center Hull (center of volume) 

P2 ■ Reference Center Tail 

P3 - Reference Center LPU 

F4 - Attach point of LPU 

P5 - Attach point of Landing Gear 

Pg - Center of gravity hull 

P7 ** Accelerometer location 

Pg - Velocity sensor location 

P9 - Rotor hub 
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Figure 2a. Fa&lc LPU Geometry Configuration Vectors 




Section d-d' 


Right Side View 


Figure 2b. Basic Landing Gear Configuration 


RA0P2 


TR-1151-2-III 


9 


PjQ - Propeller hub 

P|X ~ I^U center of gravity 

P^2 “ LPU aerodyoaalc center 


ORIGKML PAGE 
OF POOR QUALITY 


Figure 3 Illustrates the angular orientation of the LPO reference 
axes with respect to the hull. Figure. 4 shows the propeller and direct 
thrust vector orientation with respoct to the LPU. 

This progrM was designed to eodel a heavy-llft alrsulp with four 
lift propulsion tc'lts. By choosing special geoa^try and eass character- 
istics, the slaulatJ.on will aodel fewer than four LPUs. [It If sot 



Figure 3* LPU Euler Angles, 

Namelist NGBANG, Data File GMDTA, Page Volume IV 
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T Propeller 



Figure 4a. Propeller Shaft Orientation Angles (Naaellst MPRPItlG), 
Data File 6MDTA. Page A-6. Voluae IV 


T«| Direct 



Figure 4b. Direct Third Orientation Angles (Namelist NJETHSA), 
Data File GNDTA, Page A-11, Volume IV 
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pMSlbl* to Bod*! aor* than four LPUa^] This i« doa« by tarolng th« 
••rodynaalc affects of one or aore *>t the lift propulalon units and 
placing the LPUs at the vehicle center of gravity as follows: 

a) Set the nuaber of ..PU ^ to the correct nuaber 
[variable NUMLPU]. 

b) Position the LPU(s} la seroed so that their 
centers of gravity lie oci the center of gravity 
of the hull* e.g., to elialnate LPU-1, 


-RHULCG •»> RATCHl - RLTCTl t- RGGLPl 


0 

0 

0 


will place the center of .gravity of LPU-1 on the 
center of gravity of the hull* 

c) The aasa of the LPUs can never be set to sero* 
Consequently* the user should calculate the aass 
of the hull so that the aasses of the LPU(s) 
being seroed plus the Basii of the hull equals the 
desired aass of the huJ 1* The aoaents of Inertia 
aust be adjusted so the-' add in a siallar aanner. 

d) There are several aerodynaalc input*', contained 
in data file ARODTA* which aust a*.^o be set to 
sero. 

The aoorlng point geoaetry inputs are not used in the unaoored 
programs (HLASIN and HLAPAY) but they aust be in the input list* The 
landing gear spring constants and attach points are only used by HLASIM 
or HLAPAY if Che vehicle touches the ground durlni :.ne tiae history (the 
vehicle can be flown down to ground contact during a time history)* 
They are* of course* a primary factor in HLA ,R* 

The program was developed with four 'anding gears. As with the lift 
propulsion units* one or more of those landing gears can be removed* 
giving Che oftect of fewer Chan f<r . [no greater than four landing gears 
can be modeled]. To remove a landing gear* no special geometry config- 
uration is necessary; only the gear spriny constant [GEARK] should be 
set to sero for each landing gear that 1) being seroed. (l.e.* GEARKl ■ 
0.0 will eliminate all forces s'isocl .ted with ^war 1). 
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1. Dm FiLB umt iMoatk 


Contents: The Aerodyoanlc Constants for the Hull end Lift Propulsion 
Units 

These input variables can be found in th^ output listing under the 
sections headed **LPU AERODYNAMIC PARAMETERS INPUT*' and "HULL AERODYNAMIC 
PARAMETERS INPUT." 

The user will have to consult the Technical Manual to find the 
■ethod for estlnating these aerodynanlc paraaeters. Appendix A contains 
a list of all of these paraaeters and their definitions as well as a 
default value which will effectively cancel this paraaeter froa the cal- 
culations of the prograa. This default value can be used if the user: 
1) does not think that this particular effect is laiportant to be aod- 
eled; or 2) is unsure of the aethod for generating a correct value. All 
Inputs are in naaelist foraat. 

The discussion of data file GHDTA contains a dlagraa of the geometry 
configuration vectors. This dlagraa Includes the vectors, RACLPl-4, 
which are in file ARODTA. 

When the user is simulating a configuration that has less th<^n four 
lift propulsion units there are several variables in this file which 
must be zeroed. They are: the lift curve slope, the drag equation 
coefficients, and the LPU fuselage aerodynamic Y-force, Y-force, and 
Z-force derivatives. [The user should see the discussion of the geome- 
try input file, 6MDTA, for a discussion of the geometry configuration 
necessary to go with these aerodynamic parameters when zeroing some of 
the LPUs.] For example, the aerodynaaic effects of LPU-1 will be zeroed 
by setting the following values to zero: 


DLTPIA 

m 

0 

DLTRIA - 

0 

DLTPIB 

- 

0 

DLTRIB « 

0 

DLTPIC 

- 

0 

DLTRIC - 

0 

LCSPl 

m 

0 

LCSRl - 

0 

XUUAFl 

m 

0 



YVVAFl 

m 

0 



ZHWAFl 

m 

0 
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C. DATA FILS NAMB: IFCDIA 


Contents: Aerodynamic Interference Constants 

These Input variables can be found on the output listing under the 
titles: 

INTERFERENCE CONSTANTS ON ROTOR 
INTERFERENCE CONSTANTS (HI PROPELLER 
INTERFERENCE CONSTANTS ON FUSELA(X 
INTERFERENCE CONSTANTS (HI HULL 
INTERFERENCE CONSTANTS ON TAIL 

The user will have to consult the technical manual to see how 
these constants are estimated. Appendix A lists these variables with 
notes on valid input and input default values. All inputs are in name** 
list format. 
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D. Dm FILE NAME: PLMDI4 


Contents: The Rotor and Propeller Spin Rates and the Mechanical Con- 

trol Llolts on the Rotors, Propellers, and Tall 

These variables are on the output listing under "ROTOR AND PROPELLER 
SPIN RATES" and "MECHANICAL FLIGHT CONTROL SYSTEM CONSTANTS." 

The mechanical flight control system limits represent physical 
limits for the movement of the various conLrol element effectors. Data 
file HISDTA contains control system limits which correspond to these. 
These limits restrict the effector commands generated by the flight con- 
trol system loops. The flight control system loop limits should. In 
general, be less than physical (mechanical) limits, (an individual 
flight control system loop should never generate a command which uses 
all of the control ppwer possible In any particular effector). 

Appendix A contains a description of these variables. They are all 
In Namelist format. 


TR-1151-2-III 


15 



B. DATA FILE HAMB: TBMDIA 


Coatents: Trim States, Ataospherlc Paransters, aod Stability 

Derivative Flags 

These Inputs are found In the outut listing under 
e "INERTIAL VEHICLE STATE INPUTS," 
e "ATMOSPHERIC PARAMETER INPUTS," 
e "STABILITY DERIVATIVE FLAGS." 

All Inputs are In naoiellst format. Appendix A contains a description of 
each of these variables. 

This data file contains the basic data defining a trim condition. 
The number of trim runs Is defined by the response to the question "How 
many trim flight -conditions?" (see Section 2-1). The program Is 
designed so that It will read one block of trim data, do the associated 
trim calculation and print the results. It will then read a second 
block of trim data, and a third, and a fourth, up to the number of trim 
conditions requested. The data set In this file must be repeated enough 
times to prevent the program from trying to read beyond the end of the 
data file. If this data file contains more blocks of data than the num- 
ber of user requested trim conditions, the program will ignore the extra 
blocks. The last condition calculated Is used as the basic starting 
condition for the time history run which Is to follow. 

Note: 

For HLASIM and HLAPAY, the hull position (HULPOS) 
must be such that the landing gears, gear frames, 
hull belly, and tall are not touching the ground . 

This data file Is used by all three programs (HLASIM, HLAPAY, and 
HLAMOR). To Insure that the file Is correctly positioned, there are two 
flags which the input routines look for: 

1) SUBROUTINE IN3XAT must Immediately precede 

SNINSIaT each time this namelist appears. 

2) SUBROUTINE INATMOS must Immediately precede 

$NINATMOS each time this namelist appears. 
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Following each trla, the user has the option of having open "■loop 
etability derivatives calculated* and aay choose which stability deriva- 
tive oatrices are to be printed* These calculations are unaffected by 
the closed-loop control systea* The reasons for providing these options 
are twofold: 


a) The cost of calculating stability derivatives can 
be quite high because the basic program state 
derivative vector is calculated three times for 
each stability derivative matrix column* 

b) The listing of stability derivatives is rather 
long* If they are not wanted the user will pro- 
bably not want to have them printed* 

The meanings of the stability derivative flags are: 


DERVFL - T 

• F 


AMATFL - T 


BMATFL - T 


BPMTFL - T 


CMATFL - T 


CFMTFL ■ T 


Stability derivatives will be cal- 
culated* 

No stability derivatives are calcu- 
lated. [The remaining flags have 
no meaning to the program.] 

The A matrix will be calculated* 
The A matrix is generated by per- 
turbation of the state vector* 

The B matrix will be calculated* 
The B matrix is generated by per- 
turbation of all the Individual 
control element effectors* 

The B' matrix will be calculated* 
The B' matrix is generated by per- 
turbation of the six linked con- 
trols* 

The C matrix will be calculated* 
The C matrix is generated by per- 
turbation of all n'ae gust states* 

Will cause the auxiliary matrix 
associate! with each of the above 
matrices to be printed. 


Associated with each matrix there is an auxiliary matrix, which 
Is produced by the constraint forces and moments. It Is calculated slm- 
’’Itaucously with the main matrix* If the user Inputs CFMTFL ■ T, then 
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for each of the aain aatrlces* which he requested » the auxiliary aatrlx 
will also be printed. In the following exai^ile the user will have the A 
■atrix, the C natrlx» the A auxiliary aatrlx, and the C auxiliary aatrlx 
printed: 

DERVFL - T 
BMATFL - F 
BPMTFL - F 
CMATFL • T 
CFMTFL ■ T 

If CFMTFL had been false only the ' A and C aatrlces would have been 
printed. Whether the user sets CFMTFL equal to true or false has no 
effect on the cost of running the program, since auxiliary aatrlces are 
a byproduct of the main matrix calculation. 
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r. MTA PIU NAMB: HISOXA 


Concents: All Inputs Related to the Time History Calculations 

These Inputs can be found In the output listing under the headings: 

• "FLIGHT CONTROL SYSTEM CONFIGURATION" 

• "TIME HISTORY PROFILE" 

• "GUST INPUTS" 

• "COMPUTER ALGORITHM TIME STEP INPUTS" 

This file Is read only If the user has requested a tlae history calcula- 
tion In response to the Interactive, question "Tine History?" If the 

user does not want a time hlstocy calculation. It Is not necessary to 
load this file with the program. All Inputs are In namelist format. 
Appendix A contains a description of each of these variables. 

Generally, no single loop of the flight control system should com- 
mand the maximum on any Individual control effector, and no Individual 
integrator should command Che maximum cf any particular control loop. 
The Integration limits (UILM, VILM, HDTILM, PHIILM, THEILM, and RILM) 
should be less than Che corresponding loop limits (ULLM, VLLM, HDTLLM, 
PHILLM, THELLM, and RLLM). Both these sets of limits should be smaller 
than the mechanical flight control limits Input In data file PLMDTA. 
These mechanical limits represent maximum physical limits. 

The control system feedbacks can be based on two separate reference 
axes. The feedback control flags (UFDBK, VFOBK, and RFDBK) when set to 
true cause the feedbacks In the u, v, and r loops Co be based on body 
axis reference. When set to false, the u and v feedbacks are based on 
the body axis apparent velocity [relative wind], and r Is the Euler yaw 
rate, 

Any of the six control loops may be open (Inactive) or closed 
(active). The closed-loop flags are ULPFLG, VLPFLG, HOTFLG, PLPFLG, 
QLPFLG, and TRTFLG. [T Indicates Che system Is closed loop; F Indi- 
cates the system Is open loop.] The user must refer to the technical 
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manual for a description of the flight control system gains and how they 
are calculated. 

The flight control system has a hover control command. This command 
consists of the starting and ending times [POSHTl, POSHT2] for the vehi- 
cle to hold a hover position. The user does not define the position, 
but only the times. The hover position Is flagged as the Inertial posi- 
tion of the vehicle at the time POSHTl. The hover control functions as 
follows: 

1) For POSHTl < TIME < P0SHT2 the hover control 
module generates commands to hold and/or return 
the vehicle to the position flagged (PHRF:X, 

PHRF:Y, PHRF:Z, and PHRFtPSI on the output list- 
ing). These commands will react to any vehicle 
motion away from this reference position. All 
user Input flight control commands are Ignored 
during this time span. 

2) "TEST" commands will be recognized while the 
hover control Is turned on. As soon as the vehi- 
cle moves away from the hover position as the 
result of the "TEST" command the hover control 
will react to counter that motion. It will even- 
tually overcome the effect of the "TEST" command. 

3) When the hover control Is turned off (P0SHT2), 
the system resumes the input conm.and Interpola- 
tions. If the user has Input a command at the 
same time the hover control Is turned off, the 
system Immediately jumps to that command, creat- 
ing a step command. If there are no commands 
Input at the same time the hcver control Is 
turned off, the system will Interpolate between 
the last hover control command (Internally gener- 
ated) and the next flight control command (user 
Input ) . 

The user should refer to the flight control system commands below 
for further discussion and diagrams on how these two command systems 
Interact. 

The TIME HISTORY PROFILE contains the test commands and flight con- 
trol system commands which will be executed during the time history run. 
The Test Commands control the Individual element effectors; they move 
the effector a specified angle for a given Lime. [Note: T1 and T2 In 
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the variable name indicate starting and ending time, respectively*] 
Each group of effectors (l*e., rotors, propellers, and tail) have a 
single user selected starting and ending time associated with Its com- 
mand Inputs. The linked controls likewise each have one starting and 
ending time. Each of the six linked controls command one of the six 
degrees of freedom of the vehicle. These controls will manipulate all 
the various effectors according to the odxlng scheme which has beeu 
established Internal to the program. If the user wishes to change this 
mixing ocheme, the program code must be altered (see Programmer's Han- 
ual/ This should be done with great care because It may very adversely 
affect the trim and closed loop response. 

The test commands are added to the flight control system commands. 
If the vehicle Is being flotra closed loop, these test commands serve as 
control system excitations. 

The FLIGHT CONTROL SYSTEM COMMANDS are Input In a format different 
than the vither Inputs In this data file. This Is done to allow a vari- 
able number of commands to be Input without using counters or flags. 
The user may Input from zero to twenty commands for each of the six 
loops [linked controls]. Each command Is a time-command pair (real num- 
bers). It Is recommended that the commands be put on separate lines, 
although 1c Is not necessary. The time Is first and the command Is 
second. E'ch number must be followed by a comma. 

The six command variable names. In order, are: UCMD, VCMO, HDTCMD, 
PHICMD, THECMD, and TRTCMD. The namelist, NCOMAND, Is followed by 
groups of commands eacn preceded by Its variable name. If there are no 
commands for a particular variable, then that variable should not 
appear. If the user Is not Input! ng any commands for the time history, 
the namelist name, followed Immediately by the end flag must still be In 
the file. 

l.e., l$NC.iiMANO, 

0$END 
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A typtcd' ««t of coasaDd* could b«: 

l$NrOKiAND 

UCtO • 0.0, 2.0, 

3.0, 8.0, 

PHICK? > 3.0, 0.02, 

TRTCMD - 0,0, 0.0, 

2.0, O.OS, 

0$END 

Note: 

1) All of the Input veluec are followed by coaaaa. 

2) Variable naaes of controls not being cooBModed 
are excluded [l.e., VCHD, HDTCMD, and m ^]. 

3) There can be up to twenty tlM*>coiaiMnd pairs for 
each control. 

4) The tla.is [first nuaber of each pair] oust be In 
ascending order for each control. No two times 
for the ssM control can be equal. 

5) The cosmand pairs are written on separa*. a rows. 

This Is not necessary, but Is highly recoaaended. 

If the user does not Input a flight control system command at time 
0.0 for any of the six commands, the program generates a tl^e 0.0 com- 
mand equal to the trim condition. After the time history has begun, the 
program repeatedly interpolates between tho Input coasaands (or between 
trim and the first user command) to get an appropriate command for any 
particular time. In the example above, the user has Input a command of 
2*0 at time 0.0 for UCND. The program will calculate a coaunand of 4.0 
at time 1*0 [Interpolation between time 0.0 and time 3*0]* If the user 
has Input a command at time 0*0 different than the trim condition, the 
user input command is Implemnted at time 0*0* This results In a step 
command at time 0.0* 

If the user commands do not ''ontinue until the time history Is com- 
plete, the program will hold the last user command until the end of the 
simulation* This Is done to be consistent with the gust Inputs (see 
below). If the us*r commands continue beyond the final simulation time, 
the progr.''m will Interpolate toward those commands* 
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ORIGINAL PAQ:^ [3 
OF POOR QUALITY 

As noted earlier, the user can create a step conmand at tine 0.0 by 
inputting a command which is different than the trim value. A step com- 
mand can also be created during a tiaie history run by inputting two com- 
mands with times which are very close together (the times mu^t never be 
equal). For example: 

UCMD - 3.0, 25.0, 

« 3.001, 35.0, 

This commands 25 at time 3.0. The next algorithm step the coousand will 
be 35 which is, in effect, a step command. 

The following example tolll demonstrate these options: 

VHUL - 18.0, 0, -1.0, (from file TRMDTA) 

HULEUL » 0.0, 0.1, 0.0 

1§NC0MAND (from file HISDTA) 

UCMD - 3.0, 20.0, 

6.0, 30.0, 

8.0, 50.0, 

THECMD - 0.0, 0.0, 

5.0, 0.0, 

7.0, 0.15, 

12.0, 0.5, 

$END 

TSIM - 10.0, (simulation time; file HISDTA) 

Note: -1.0 in the trim velocity, VHUL * 18.0, 0.0, -1.0, refers to 

a flight condition of 13 units forward (positive) and 1 unit vertical 
(upward) velocity. This negative body axis velocity corresponds to a 
positive climb rate command (HOTCMD > -t-l.O). 

Figure 5 Illustrates the example above. Figures 6 and 7 illustrate 
the operation of the hover control system. 

" le GUST INPUTS are next in this data file. There are two types of 
gust inputs, a (l-minus-cosine) gust generated at individual vehicle 
elements, and four gust strings which are read from four user created 
files. 

The beginning of these gust Inputs must be flagged with SUBROUTINE 
INGUST immediately preceding namelist NHGCOM. This is done to allow 
HLAMOR to skip over the unused flight control system and time history 
profile Inputs. 
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ORIGINAL PAGC 13 
OP POOR QUAUTV 


X Trim condition 
O User input commond 
• Progrom generated commands 
— — Commands generated by interpolation 


UCMO 

(ft/sec) 



3 - 

HDTCMD 2 - 

(ft/sec) I ^ 

0 I 1 1 1 1 1 L 

0 2 4 6 8 10 !2 

Time (sec) 


THECMD 

(rad) 



Figure 5. Program Cenerated Control Commands Eased 

on User Input Commands 
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^^ginal pass ra 
POOR QUALITY 


Flight Control Contmonds input 
UCMO • 0.0 , 0.0, 
2 . 0 , 00 , 

6.0 . 15.0, 

10 . 0 . 20 . 0 , 


Hovor Position Commond 
POSHTI ■ 2.0, 
P0SHT2 * 6.0 



Time (sec) 


O User input command 

System calculated commands by interpolation 
between user inputs 

Commands which would hove been used If hover 

control had not been turned on 
/'\y\ Commands which the hover control generates 


Ftjiuro b. Hover Control Operation-Hovering 
Trim Condition 




2 ') 



Commands (ft/sec) 


ORIGINAL PAGE fS 
OF POOR QUALITY 


Control C'>mmQnds input 
UCMO « 0.0, 10.0, 
2 . 0 , 10 . 0 , 

6 0 , 20 . 0 , 

10.0, 15.0, 


Hover Position Commond 
POSHTI • 3.0 
P0SHT2 « 5.0 



Figure 7. Hover Control Operation- 
Nonhovering Trim Condition 
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Tht (1 ’■minus -cosine) gust on eech element Is turned on end off by 
the starting (Tl) and ending (T2) time, similar to the format of tha 
test commands and the hover control command. These gusts trace out a 
(l-mlnus-cojlne) curve over the time span of the gust. This is dons by 
Inputting the gust time span (variables with Tl and T2 in them) and the 
maximum gust values (variables* with MAX In them). If a (one-minus- 
cosine) Input Is not wanted on any or all of the vehicle components* the 
' ' irtlng times should be set to values larger than the simulation time* 
or the maximum values of the gust velocities should be sero. The ending 
time for the gust must always be larger than the starting time. There 
are 6 points on the vehicle where these gusts can be generated: Hull* 

Tall* and each of the four LPUs (see Figure 8). 

The starting and ending times for the various gusts can be differ- 
ent. This allows the user to stagger the starting and ending times in 
such a manner as to approximate a guec field passing by the vehicle. 
This is only a crude approximation and is not really the intention of 
the model. The Intention Is to allow the user to Input a gust which 
acts at only one point on the vehicle and therefore isolute the effects 
of that gust. 

The gust string Input option Is turned on by the flag GSTFLG > T. 
The user must load the gust string Input files (RG1-RG6) with the pro- 
gram only if this flag Is turned on. Discussion of files (RG1-RG4) is 
found In Section 2-G. RG5 and RG6 are for payload input files gusts and 
are discussed In Sectli>n 7* Subsection B* Article 4. The three vari- 
ables, RFSRCX, RASRCX and RSORCX* determine the position of the gust 
input sources around the vehicle (In vehicle reference frame). The two 
x-varlables* R'p'SRCX and RASRCX, define the forward and aft vectors to 
the Input source. The single variable RSORCY (always positive) Is the 
distance right and left to the input source locations from the hull cen- 
ter line. This insures that the four gust Input sources are located in 
a laterally symmetric rectangular configuration about the vehicle. It 
Is recommended that the rectangle defined by those sources should 
surround the LPUs and tall, but not be unnecessarily large. Figure 9 
Illustrates a typical positioning of the sources about the vehicle. The 
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(I /sec) (ft/sec) 


ORIGINAL PAGE « 
OF POOR QUALITY 


$OT1GC0M 

HTIGST 

HT2GST 

TJHCMAX 

VHGMAX 

WIIG'IAX 

PHGMAX 


3.0 seconds 

7.0 seconds 
0.0 ft/sec 

8.0 ft/sec 

3.0 ft/sec 
0.0 rad/sec 


QHGM/OC 

RHGMAX 

nUXHMX 

DIJYHMX 

OVYPKX 

SEND 






dv 

^ Gust Velocity Gradient 
dy 


Figure 8. Typical One-Minus-Cosine 


0.0 rad/ sec 
0.0 rad/ sec 
0.0 1/sec 
0.0 I/sec 
10.0 1/sec 


Gust z- Direction 



Time (sec) 


Gust 
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Source I 
Source Z 
Source 3 
Source 4 


^fGlNAL 
Of POOR 


iil 

quality 


Source I 


Source 2 


RGI 

R62 

RG3 

RG4 


Source 3 



RFSRCX 


RASRCX 


Figure 9. Gust Input Source Geometry; RFSRCX Is 
Positive as shown. RASRCX is Negative as 
shown. RSORCY is Always Positive 
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lnterpol«tlon «ch«M trtilch u««« Ch««« aourc«« do«« not conaidar tha dll'- 
farant halghta of tha varloua vahlcla coaponanta* 

Tha data undar “COMPUTER ALGORITHM TIME STEP INPUTS" ara tha Uat 
Inputa froa thla flla. Thay ara contalnad In Naaallat NINSTBP, which 
auat ba pracaded by tha flag SUBROUTINE INSTEP. 

TIHSTP la tha tiaa atap that la paaaad to tha ayataa Intagrator* 
Tha Intagrator dlvldaa thla valua Into aaallar and aaallar dlvlalona 
until it can aaat tha Intarnal arror criteria. Thla variable provldaa 
tha Initial atapalia to tna Intagrator. It la alao tha naxlaun atap tha 
Intagrator can taka. TIHSTP ahould ba raaaonabla for tha vahlcla aoda 
fraquanclaa; unnacaaaarlly aaall valuaa will Incraaaa program axacutlon 
tlaa and coat. During davalopment wa found that valuaa batwaan 0.1 and 
O.S tiorkad wall. 

Tha variable. MINSTP la tha olnlaum allowable tlam atap which tha 
aystam intagrator can taka. If tha ayataa intagrator la unable to meat 
tha arror criteria without going below this tlaa step, the values of tha 
last actaapt ara accepted and tha program calculations continva after a 
warning message has bean printed. 

Factors affecting the choice of TINSTP and MINSTP: 

1) Rapidly changing flight conditions require saal'- 
ler minimum intagrator time steps (small MINSTP). 
l.a.. gusts, and flight control comaands. 

2) Springs In payload cables and landing gears 
require smaller minimum steps (MIKSTF' because of 
high frequency modes. (HLAPAY needs smallest 
steps. HLASIM Is only affected If active landing 
gears touch the ground). 

3) If the early part of the tloM history has chang- 
ing flight conditions but the remainder la In 
steady flight MINSTP should be small, but TINSTP 
considerably larger to allow the Integrator to 
open up Its calculation step. 

A) The plotting files (see Section 3-B) are written 
every TIMSTP interval of the time history. 
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G. PILB HAMBSs IGl, 162, SC3, IGA ^ 

ORIGINAL R;g£ ; i 

Contents: The gust Input stilngs ^ POOR O.iALJTY 

The gust input strings are read from these four data files. These 
files contain four numbers on each row. The first number is time, fol- 
lowed by the u, v, and w (linear) velocities of the gust. These numbers 
are input in free format with an unlimited number of lines. The only 
restriction is that the times (first entry on each line) be in ascending 
order. Note that first gust time entry need not be 0.0. 

Gust velocities at the source location are obtained from time inter- 
polation between the times in these files. The gust velocities and gra- 
dients for each element are determined from spatial Interpolation among 
the pertinent source velocities. Backward difference schemes are used 
to determine gust accelerations. 

These four data files are associated with the four gust input sour- 
ces (see data file HISOTA). A typical set of gust input data and the 
associated explanation is presented in the following example. 

Sample data files: 


RGl 




RG2 




0.0, 

5.0, 

0.0, 

0.0 

0.0, 

0.0, 

2.0, 

0.0 

1.0, 

10,0, 

0.0, 

0.0 

l.C, 

8.0, 

5.0, 

0.0 

^.U, 

8.0, 

5.0, 

-2.0 

1.5, 

10,0, 

6.0, 

0.0 

3.0, 

2.0, 

1.0, 

0.0 

2.0, 

12.0, 

3.0, 

0.0 





3.0, 

8.0, 

0.0, 

0.0 





4.0, 

0.0, 

0.0, 

0.0 

RG3 




RG4 




3.0, 

0.0, 

0.0, 

0.0 

3.0, 

1.0, 

0.0, 

0.0 

4.0, 

4.0, 

0.0, 

0.0 

3.5, 

6.5, 

-2.0, 

0.0 

5.0, 

0.0, 

0.0, 

0.0 

4.0, 

8.0, 

-5.0, 

0.0 





S.O, 

15.0, 

-6.0, 

0.0 





6.0, 

18.0, 

-8.0, 

0.0 





7.0, 

12.0, 

-2.0, 

0.0 


Notes: 

1) the various times (first entry of each line) on 
the different files need not correspond with each 
other. 
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2) In file RG1« there will be a u-gust of 5.0 at 
tine 0.0. This gust will not appear on the trla 
data or affect the trla/stablllty derivative cal- 
culations. It will appear at tiaa 0.0 on the 
printout and will be used in the time history 
solution for time 0.0. 

3) The user» if he wishes to have the gusts set to 
sero after a certain point* aust input sero gust 
values at that time. The prograa does not sero 
them off at the end of the data file* but rather 
continues the last input indefinitely. For exam- 
ple * if these gust files were being used for a 
10 second tiae history* source RGl %fould continue 
to have a gust velocity vector of (2.0* 1.0, 0.0) 
from the third second through the end of the tiae 
history. Source 2 (RG2) will have (0.0* 0.0* 
0.0) gust velocity after time 4.0 seconds. 

4) In files RG3 and RG4* the first line of input has 
a time of 3.0. The prograa* will autoaatically 
set the gust velocities to sero at tiae 0.0 and 
will then proceed to interpolate between sero and 
the first user input gust velocity at* in this 
case* three seconds. The velocity vector of 
(0*0*0) at 3.0 seconds in file RG3 causes no gust 
to appear at source 3 until after the third 
second. In file RG4* there is a gust of 1.0 at 
three seconds causing gusts to appear launedlately 
after Time 0.0 [l.e.* 1.0 second has (0.333* 0.0, 
0.0) and 2.0 seconds has (0.667, 0.0* 0.0)]. 
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ORiGJNAL PAGE B 
OF POOR QUAUTV 


H. DATA FILE: ODTLST 

Contents: The code numbers Indicating which variables the 

user wishes to have printed on the output listing 

This flla consists of two lists of Integers. The first list con- 
ta_ 18 the code numbers for the hull output variables, and the second 
contains the code numbers for the lift propulsion unit output variables. 
Each list must be terminated by a zero or negative number. Appendix 0 
contains a list of these variables and their code numbers. The order In 
which the code numbers appear In the Input file determines the order In 
which the variables will be printed on the output listing. A variable 
can be requested twice and It will be printed twice In the order In 
which It Is listed. The following example shows the Input and the asso- 
ciated output. 


file 

001 

OUTLST 


Resulting Output 



002 


u 

V 

w 

p 

Q 

R 

003 

004 

HULL 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

005 


PHID 

THETD 

PSID 

u 

V 

w 

006 

LPUl 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

000 

LPU2 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

004 

LPU3 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

005 

006 
001 
002 
003 

LPU4 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 
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I. FILE NAME: INPUT (TAPES IMTKBACTIVE IBSPOHSBS) 


The program waa developed on an Interactive eyatem with five quea- 
tlona to establish the basic run parameters* The user must input the 
respond - to these five questions in the data file INPUT (TAPES) and load 
it with the program. The questions are: 

1. "Six degree of freedom simulation? T/F" 

T ~ The program will continue after completing 
the trim and stability derivatives to run a time 
nlstory. The time history uses the data read 
from data file HISDTA, aud will start with the 
last trim condition calculated. 

F ~ The program will stop after completing its 
trim calculations. File HISOTA is not necessary 
in this situation* 

2. "How many trim flight conditions?" I^e answer 

must be a positive non zero Integer in the first 
two columns (12 format). This number must be 
less or equal to the number of sets of trim data 
found in file TRMOTA. If the number of trim data 
sets is more than this number, the extra trim 

sets will be Ignored. If there are not enough 

trim sets to match this number, the program will 
abort for attempting to read beyond the end of 

file. 

3. "Generation of plotting files? T/F" 

T — All output variables which are calculated 
during the program run will be written to a file 
(PLOT) at the end of each algorithm step. This 
file is to be read by a post processor program 
called PPLOTF. 

F ~ No plotting files will be written. 

4. "Do you want English units? T/F" 


T — English units will 

be printed with the 

input 

variables. 

horsepower. 

The 

power 

requirements 

will 

be 

In 

F ~ Metric 

units 

will 

be printed with the 

Input 

variables. 

Kilowatts. 

The 

power 

requirements 

will 

be 

In 
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NOTE; No aunerical convrsioas «re done to the 
Input data. Thl flag only cauaea the prograa to 
print either aetrlc or English unite, and calcu> 
late the power appropriately. 

5. "Pull header T/F" 

T — > The output listing will Include a complete 
Hat of the Input parameters. 

F -** Only an abbreviated form of the output list** 
Ing will be printed. Host of the input variables 
are not on •’*'18 listing. All output variables 
and calculated data Is Included In this abbrevi- 
ated listing. 

6. "Any consents? (6 lines)" At this point the user 
may input six lines of alphanumeric Information. 
This Information will be printed at the top of 
the output listing and Is Intended to allow the 
user to record some descriptive comments about 
the particular simulation run. 


TR-1151-2-III 


35 



J. DATA PILE: BIMSS6 


ConCtnttt Error McttagM 

This data fila la auppUad with tha prograar It la not producad by 
tha uaar. It contnlna all arror aaaaagaa which nay ba prlntad by any of 
tha prograaa* Bach naaaaga la threa llnaat 

1) Coda nuobar (i3 Format ) 

2} Up to 80 charactara of tha aaaaaga (8A10 Format) 

3} Laat 40 charactara of tha aaaaaga or a blank llna 
If not naadad (4A10 Format)* 

Tha uaar ahould la'ar to tha Ptoaraaaara* Hanual Error Procaaalng 
section If it ia nacasaary to axtand this fila* Appendix B of this man- 
ual contains a list of tha arror massagaa and a brlaf axplanatlon of 
each. 

Data file ERMSSG follows: 


001 

ABSOLUTE VALUE OF PANGLE IS GREATER TOAN 1/2 PI. 

002 

CONTROL COMMANU TIMES WERE NOT INPUT IN INCREASING ORDER. 

003 

CONVERGED SOLUTION OF CT AND WIN 13 INCORRECT. 

004 

CT AND WIN DID NOT CONVERGE. 

005 

TVe COLUMN NUMBER EXCEEDS MAXIMUM ALLOWED. 

006 

rve ROW NUMBER EXCEEDS MAXIMUM ALLOWED. 

007 

SRCWN WILL EXCEED 30. 

008 

GUSTTl IS GREATER-EQUAL TO GUSTT2. 
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009 

LCS OR SIGMA IS LESS THAN ZERO 


^IGINAL page 1$ 

OF POOR QUALITY 


010 

LENGTH OF VCTR IS NOT 6, 12, 24 OR 42. 

011 

MORE THAN 20 CONTROL COMMANDS WERE INPUT. 

012 

NO REAL POSITIVE ROOTS HERE FOUND BY THE IMSL BOUTINE. 

Jl3 

CURRENT AERODYNAMIC ANGLES DO NOT SATISFY ANY OF THE POSSIBLE CCHIDI- 
TIONS. 

014 

SQROOT IS NEAR ZERO. POSSIBLE DIVISION BY ZERO. 

015 

STABUITY DERIVATIVE HILL NOT BE CALCULATED FOR THIS TRIM. 

016 

STALL REGION ANGLE 1 IS GREATER THAN STALL REGION ANGLE 2. 

017 

SOME OF THE STALL REGION ANGLES WERE NEGATIVE 

018 

SOME OF THE AERODYNAMIC ANGLES OF THE TAIL ARE GREATER THAN PI. 

019 

TIME IS GREATER THAN LAST COMMAND TIME HHI^H SHOULD BE THE SAME AS THE 
FINAL SIMULATION TIME. 

020 

TIME IS LESS THAN THE FIRST COMMAND TIME WHICH SHOULD BE ZERO. 

021 

TICOM IS GREATER-EQUAL TO T2C0M. 

022 

INCORRECT INPUTS 

023 

IMSL ROUTINE HAS RETURNED AN ERROR FLAG. ROUTINE NAME IS THE FIRST 
VARIABLE GIVE*' BELOW. 

024 

CDFLAG IS NOT SET TO -1, 0 OR ’ ON RETURN FROM SUBR. ITERCT. 
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025 

LESS THAS 4 ZEROS WERE FOUND BY DfSL ZRPOLT 

026 

REQUIRED TRIM CONTROL EXCEEDED AVAILABLE INTEGRATION LIMITS. IF LOOP 
CLOSED THE INTE<a(AT(»l HILL BE SET TO LIMIT 

027 

THE TINE HAS MOVED BACK FURTHER THAN THE RANDOM GUST BUFFER HAS STORED 
GUST VALUES. REDUCE TINSTP ON INKJT. 

028 

THE TIME READ FROM THIS FILE IS LESS-EQUAL TO ZERO. THE TIME AND GUST 
VELOCITY WILL BE IGNORED. 

029 

THE TINE IS OlEATER-EQUAL TO 100000. 

030 

CONDITION FLAGS FROM IMSL ROUTINE DVERK. 

031 

TIME INCREMENT IS LESS THAN ZERO 

032 

THE LENGTH OF THE VECTOR PASSED INTO PPTURB IS NOT 6 (» 12. 

033 

THE VALUE OF VCTRFL IS NOT VALID 

034 

SOME OF THE INVALID STABILITY KRIVATIVES HAVE NOT BEEN FLAGGED BECAUSE 
THE ARRAY IS FULL. 

035 

THE LINEARIZATION LINEAR INCREMENTS ARE LARGE ENOUGH TO CAUSE SOME OF 
THE CABLES TO GO SLACK. THEY ARE BEING RESET 

037 

THE LENGTH OF THE SV VECTOR IS NOT CONSISTANT WITH THE SIZE OF THE BLANK 
BLOCKS FOR EXTRA INTEGRATOR STATES. 

038 

THE TIMSTP OR MINSTP INPUT IS GREATER THAN THE APPROX. PERIOD/IO AND MAY 
CAUSE NUMERICAL INACCURACIES. 

039 

LMSL DVERK WAS UNABLE TO REACH THE SPECIFIED ERROR CRITERIA WITHOUT 
GOING BELOW THE MINIMUM TIME STEP. 

040 

THE FUG FOR THIS SUBROUTINE WAS NOT FOUND IN THE DATA FILE 

041 

HAKE ANGLE 1 MUST BE LESS THAN ANGLE 2, AND BOTH MUST BE BETWEEN 0 AND 
2*PI 
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042 

SOME OF THESE VALUES HILL CAUSE DIVISION BY ZERO. 

043 

HORE THAN THE MAX NUMBER OF OUTPUT VARIABLES WERE REQUESTED. 

044 

AN INITIAL GUESS WITH LANDING GEAR IN GROUND CONTACT AND PITCH ANGLE 
LESS THAN 1. COULD NOT BE FOUND. 

045 

LINEARIZATION INCREMENT GOULD LIFT ONE OF THE LANDING GEARS OFF THE 
ROUND. IT IS BEING RESET. 
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SBCTION 3 


PtOGIAM OUTPUT 


4. OUTPUT LISTUV'- 

The output listing Is wrlttsn to fils OUTPUT (TAPE6)* Ths discus- 
sion which follows assusss ths user has Input "Trus" to ths Intsractlvs 
question "Full Header?" (see ths Interactive questions section). If the 
has user Input "False" In response to this question, the output listing 
will be abbreviated. While reading this section of the oanual, the user 
will need to consult Appendix C (saeple listings) of the User*s Manual 
Appendices, (Voluae IV). 

The listing starts with a title and a block of InfortMtlon entitled 
Run Description. These six lines, describing the run which Is being 
made, are provided by the user in response to interactive question "Any 
Comaents?" This Is followed by a coaplete listing of the program Input 
variables. These values are listed with the variable naae and brief 
descriptions of each variable. A more detailed definition of these 
variables can be found in the Programmer*” Hanual Appendices. The 
sections entitled "Geometry Inputs" through "Mechanical Flight Control 
System Constants" define the basic vehicle configuration. The variables 
found In sections entitled "Inertial Vehicle State Inputs" through 
"Stability Derivative Flags" define the basic trim condition and sta- 
bility derivatives matrices to be calculated. 

Following these Inputs, the results of the trim calculations are 
printed. The "Trim Case Number" counts the trims which are calculated 
during this particular run. The next statement Indicates whether the 
trim reached a converged solution or not. If the trim did not converge, 
the data printed here represents the best solution found. The trim 
algorithm control variables may Indicate the nature of the problem If 
the trim does not converge. A complete discussion of the trim Is con- 
tained In Section 5. 
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The data from the trim calculation consists of a standard data frame 
starting with "Simulation Time", which will be a negative number* These 
negative times will count toward zero for each trim solution calculated 
(If there are 3 trim calculations, the simulation times would be -3, -2, 
and -1). The data which follows are the user selected variables (see 
Input file OUTLST). There are more than 300 variables related to the 
hull which the user may request and more than 150 variables related to 
each LPU. The user should refer to the Section 2-H for an explana- 
tion on how to choose these variables. Appendix D at the back of this 
manual contains a list of the available output variables and their code 
numbers. 

At the end of this data frame the "Tall Aerodynamic Regimes" are 
printed for angle of attack, sideslip, and rolling angle of attack. 

• "1" Indicates prestall regime, 

• "2" Indicates transition regime, 

• "3" indicates the stall regime. 

Additional messages may be printed Indicating special conditions that 
exist at this time. They Include: 

1) Vortex ring state condition on any propeller or 
rotor. 

2) Any maximum control limits which have been 
encountered (see data file PLMDTA). 

3) Any parts of the vehicle which are in ground con- 
tact; Including landing gears, landing gear 
frame, tall, belly, bow, or stern. 

4) Hover control being turned on. 

At the end of each trim data frame the stability derivatives. If 
requested, are calculated and written to the output listing. Section 6 
of this manual discusses the stability derivatives. 

If the user has input a number greater than 01 In response to the 
Interactive question "How Many Trim Flight Conditions?" the program will 
read a new set of trim input data from data file TRMDTA and calculate a 
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new trim condition. Note, the user must have sufficient data blocks in 
data file IRNHIA to match the number of trim conditions he requests* 
The output listing will repeat all the information from "Inertial Vehi- 
cle State Inputs" through the Stability Derivatives as well as the new 
trim calculation results. 

If the user has input False to the interactive question ’Time His- 
tories?" the program will now stop. If the user input True, the program 
will read the time history data from file HISDTA and write that data to 
the output listing. That data includes "Flight Control System Configu- 
ration," "Time History Profile," "Gust Inputs," and "Computer Algorithm 
Time Steps Input." 

Following the input of the above data, the program will print a data 
frame at time 0.0 (i.e., before doing any time integration). This out- 
put should be the same as the trim data frame (time -1.0), except for 
variables related to time dependent Inputs. For example, if the user 
has input a gust string with a nonzero gust at time 0.0, the gust value 
printed will not be 0.0, and the forces related to the wind will be dif- 
ferent then the trim values. Similarly, if the user has input a flight 
control cooimand which is different than the trim control at time 0.0, 
the flight control system will alter the commands, and consequently, 
there will be changes in the resulting forces at time 0.0. Test com- 
mands can also have the same effect. In all these cases the inputs can 
alter the forces but they cannot change the vehicle position or veloci- 
ties at time 0.0. They will cause the vehicle to immediately begin mov- 
ing away from the trim condition when the integration begins. 

After printing the clme 0.0 data frame, the program Integrates the 
vehicle states forward in tliK until reaching the user input print time 
step. The data frame is again printed and the process is repeated xmtil 
reaching the simulation time. The program then stops. 

If the program encounters an error condition at any point in the 
simulation, a message is printed immediately. Consequently, the error 
occurred during the following, not the preceding data frame calculation. 
If an error message appears, the user should consult the Error Process- 
ing section of this manual. 
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B. PLOTTING FILE 


If the user eaters True in response to the Interactive question 
"Plotting Files?", a binary format file (PIX)T) will be created with all 
of the output variable names and values. This data is written after 
every trim calculation and at every algorithm (TIMSTP) step (not the 
print step; see Section 2-F) during the time history run. Therefore, 
the plotting file will contain much more data than the output listing. 
All output data generated by the program (see Appendix D) are written on 
the plotting file. The code numbers in file OUTLST do not effect the 
variables written to the plotting files. Figure 10 illustrates the 
structure of the plotting files. 

The file starts with the program ID Indicating which program gener- 
ated the file. Following that, the Julian date, and then a record 
containing the number of hull variables and the number of LPU variables 
are listed. Following the above three records, the data file iray vary. 
Programs HLASIM and HLAMOR generate identical data files (left column). 
Program HLAPAY generates a file with the same structure, but intersper- 
ses the payload data with the vehicle data (right column). 

During the trim calculations, the payload data is written first. 
(Note, trim calculations are indicated by negative time). If there is 
more than one *:rim calculation, the times will count as negative numbers 
going toward zero. At the last block of trim data (time equal to -1.0) 
the variable names are printed out following the variables for both the 
payload and the vehicle. 

At this point we enter the time history section of the program. If 
the user did not run a time history, then, o£ course, the remainder of 
the data file will not exist. Immediately following the variable names, 
the program algorithm time step and the total simulation time is writ- 
ten. Then, the data for each time frame of the time history is written; 
the time is listed first followed by the variables. This is repeated 
for every algorithm step until the end of the simulation. The HLAPAY 
time history data has the vehicle data written before the payload data 
(the opposite of the trim). 
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Program ID - HIASIM, HLAHOR, or HLAPAY 
Date (Julian) 

Number of hull variables, number of LPU variables 



Program ID HLASIM or HLAMOR 

Program ID HLAPAY 



Number of Payload variables 
and number of cable variables. 



Time (Negative during Trim) 
Payload variables 
Cable variables 

Repeated for 
multiple 
trims, until 
time - -IsO 


(When Time ■ -l.O) 
Payload variable names 
Cable variable names 


Time (Negative during Trim) 

Time (Negative during Trim) 


Hull variables 

Hull variables 


LPU variables 

LPU variables 


(When Time ■ -1.0) 

(When Time - -1.0) 


Hull variable names 

Hull variable names 


LPU variable names 

LPU variable names 

Only if Time 

Algorithm time step and total 
simulation time 

Algorithm time step and total 
simulation time 

Histories 

Repeated for 
each algor- 

Time (positive during time 
history) 

Hull variables 

Time (positive during time 
history) 

Hull variables 

ithm step 

LPU variables 

LPU variables 



Time (positive during time 
history) 



Payload variables 
Cable variables 


Figure 10. Format of the Binary Plotting File (PLOT) 
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This plottiag fils is to bs ussd by s post procsssor which gsnsrstss 
filss sppropriats to ths ussr's plottiag squipownt and softwsrs. Ths 
vsrisbls nsass srs writtsn in ths ssas ordsr as ths variablss, allowing 
thsa to bs pairsd by uss of array subscripts. Appsndix D contains ths 
vsrisbls nsass and ths cods nuabsrs. By osttching vsrisbls nsass and 
vsrisblss, and using ths tias (nsgstivs for tria» positivs for tins his- 
tory) ths post procsssor can sort ths data. 

As part of ths iaplsasntation of this projsct on ths NASA/Aass 
7600/Scops Systea, a post procsssor is bsing supplisd to gsnsrats plot- 
ting filss of ths tias history data (sss Ssction 9 of this annual). Ths 
ussr aay wish to writs a post procsssor which will construct plotting 
filss for plotting ssts of tria solutions in a aannsr siailar to ths 
tias history post procsssor. 
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SBCTION 4 


HIM 


Th« trla algorlthn is a geaarallsed sacant schSM «rt)lch Iterates 
until controls are found which reduce the accelerations below the trim 
tolerance. The trimmer manipulates six general controls in searching 
for the converged solution.* They are: 


UDCNTL 

“ x-directlon control 

VDCNTL 

- y-direction control 

WDCNTL 

- z-direction control 

PCONTL 

- roll control 

QCDNTL 

- pitch control 

RCONTL 

- yaw control 

The limits 

are: 


Maximum Iterations - 200 

Maximum Restarts 6 

Maximum Trim Tolerance ■ 10 ^ (Convergence Criteria) 

The ’’Iterations" are a count of the number of trim algorithm passes 
used in converging to the solution* The "Restarts" are a count of the 
number of times the trim algorithm was restarted because it encountered 
a local solution which did not conform to convergence requirements 
(i*ee» "local minimum")* The trim is stopped when the convergence cri- 
teria is satisfied (converged), or the maximum allowable iterations or 
restarts is exceeded (i*e*, solution not converged)* The convergence 
criteria is a measure of the trimmed acceleration of the vehicle* The 
trim algorithm is based on finding controls which zero all accelerations 


^While reading the following discussions the user should refer to 
the sample output listings in Appendix C* 
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acting on the vehicle for the Initial conditions specified by the user 
on Input file TRMDTA. 

The "Trim Controls" are the six Iltdced controls irtilch achieve the 
convergence criteria for a specific flight condition* The following set 
of flags count various error conditions that nay occur during the trim 
calculations: 

XT HER - The number of times the collective pitch on 
each rc'or exceeded the user Input limits (data file 
PLMDTA) . 

TTHEP - The number of times the collective pitch or 
each propeller exceeded the user Input limits (data 
file PLMDTA). 

TAISR ~ The number of times the lateral control axis 
deflection on each rotor exceeded the user Input 
limits (data file PLMDTA). 

TBISR - The number of times the longitudinal cyclic 
pitch angle on each rotor exceeded the user Input 
limits (data file PLMDTA). 

SNGMTX - The number of times the trim alytcrlthm 

encountered singular matrices. 

AILERON - The number of times the allercn control 
exceeded the user Input limit (data file PLMDTA). 

ELEVATOR - The number of times the elevator contr ' 
exceeded the user Input limit (data file PLMDTA). 

RUDDER - The number of times the rudder control 

exceeded the user Input limit (data file PLMDTA). 

Note, the four numbers under the rotor and propeller elements refer to 
LPUs 1-4 respectively. 

These flags will sometimes Indicate the nature of the problem If the 
trim does not converge, l.e., 

THEP 0 58 0 52 

Indicates the trimmer was asking for more than the maximum collective 
pitch on LPU propellers 2 and 4 (yaw control) than available (THEPMX In 
data file PLMDTA). 
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SBCnON 5 


SZABILITT BmVATIVBS 

Following each trim solution, stability derivatives will be calcu- 
lated upon request and printed. The stability derivative output begins 
with a header which describes the basic form of the matrix equation that 
generated the stability derivatives* There are four basic matrices and 
four auxiliary matrices (the user has the option to choose any or all of 
the matrices, see data file TRMDTA). These matrices are generated by 
the forward-backward difference scheme on each of the elements of the 
vehicle state vector, the Individual control elements, the linked con- 
trols, and the gust- values at each vehicle component* The four matrices 
are: 

1) A Matrix - The A Matrix Is generated by perturb- 
ing each of the elements of the state vector* 

2) The B' Matrix - The B' Matrix Is generated by 
perturbing each of the six linked controls* 

3) The B Matrix - The B Matrix Is generated by per- 
turbing each of the control elements on the four 
LPUs (rotors and propellers) as well as the three 
nail contiols* 

4) The C Matrix - The C Matrix is ^.aerated by per- 
turbing each of the gust values on the vehicle, 
tall, and the four LPUs. 

Along with the basic stabllllty derivative matrices, a set of auxil- 
iary matrices are generated, based on the constraint forces between the 
vehicle and the four LPUs* These matrices, A Auxiliary, B' Auxiliary, B 
Auxiliary, and C Auxiliary, correspond with the A, B', B and C Matrices* 

In the p- ocess of calculating these matrices the program checks for 
stiong noullnearltles* Following the matrices, the program lists those 
elements of the matrices where nonllnearltles were detected* It prints 
the derivative at the positive Increment and the negative Increment of 
the forward-backward scheme so the user can see the nature of the non- 
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linearity* Figure 11 illustrates the three derivatives in relation to 
the forward/backward t^<’rturbations* 

The eigenvalues and normalised eigenvectors (real, imaginary parts) 
are then printed. The normalised eigenvectors are printed in coluius 
directly below the corresponding eigenvalues* The IMSL performarce 
ind is an Indication of the accuracy of theue values* A performance 
ln< ex below 1*0 indicates the accuracy is excellent* An index between 
1*0 and 100, indicates that the accuracy is good, and indices greater 
than 100 Indicate poor accuracy* The user is refered to the IMSL Manual 
(subroutine EIGRF) for a more detailed explanation of this performance 

I 

index* 

Finally, for the user's information, the perturoatlon Ir. *ements 
which were used in the calculation of the sta' ility derivatives are 
printed* These increments are set internally in the program (subroutine 
INTIAL)* 
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TR-1 

-w 


sjassw 



Backward and Forword 
Perturbation Increments 


Pq - initial values 

Pi - Result of negative perturbotion 

P 2 - Result of positive perturbation 

Slope d, - Positive Increment derivative i printed with nonlinear 

Slope da - Negative increment derivative > element 

Slope do - Derivative vaL rs in stability derivative matrix 

If |d 2 ~d|| ^0.3 do the element will be flagged as invalid 


Figure 11. Stability Derivative Calculation Algoritlun 
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SECTION 6 


ntocusiiiG 

Th«r« «r« *«v«r«l dtff«r«nt conditions which ctuss ths progrsn to 
print sn error assssge* They ere: 

1. User errors - These could be Invalid l>iputs, 
nuabers which asy cause division by teroi or 
other values which the prograa determines are 
Invalid for that particular variable* 

2. Model error conditions - Conditions with which 
rhe prograa aodels have not been designed to 
cope* They aay not necessarily be wrong, but 
only that the prograa Is unable to deal with 
thea* 

1* Conditions of particular Interest — These are 
conditions which aay be of Interest to the coapu- 
ter prograoHser or the engineer* 

These conditions aay or aay not cause the prograa to stop* If the con- 
dition does not invalidate the calculations, the prograa Is allowed to 
continue. The error aessages are Input as data Into the prograa and are 
supplied with the prograa (see Input data file ERMSSG)* Appendix E 
lists the error aessages which aay be printed with an explanation of 
what they iwan and suggestions for rectifying the problea. 

A typical error message consists of a lire of stars, followed by a 
one line error message, the subroutine name which detected the error, 
and up to three variables with their present values froa that subrou- 
tine. These variables should prove useful If It Is necessary for the 
user to consult the prograa code to ascertain the cause of the error* 
The following example shows the format of the error message. 


22 INVALID INPUTjJ 
SUBROUTINE INLAAO 
BETAIT - .6981 E+00 

BETA2T - .5236E-H)0 
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In this particular case the Input routine INLARO has detected that 
BETAlT Is not less than BETA2T as It should be* "22'* Is the error aas- 
sage code number (see data file ERMSSG). 
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SBCTIOif 7 


PROGRAM HLAPAI; SIMULATION Of A HEAVY LIFT AIRSHIP 
CARRYING A PAYLOAD ATTACHED BY GABLES 


A. INTRODOCTION 


This program models the same vehicle as program HLASIM. It has, la 
addition, the models of a slung payload and supporting (elastic) cables* 
The payload and cable models were developed on an entirely separate pro- 
gram. HLAPAY contains both the program HLASIM and the payload-cable 
program. They Interact through the cable forces. The user should keep 
this modularity In mind while creating his data file and studying the 
output from this program. The reasons for Implementing these models In 
this way are: 

1) It allows for a modular development and debugging 
of the program. 

2) Will allow future alterations of either the pay- 
load or vehicle to be Implemented and tested 
separately before being put together as one pro- 
gram. 

It Is expected that the user has already used and has a good under- 
standing of program HLASIM before attempting to use program HLAPAY. 
Consequently, In this section of the manual there %d.ll be no discussion 
of the use of program HLASIM. We present here a detailed description of 
the additional Inputs and considerations necessary to study and use pro- 
gram HLAPAY. The payload model consists of: 

1) A payload module and a cable module. 

2) Input subroutines which read and then write to 
the output listings all the payload related 
Inputs. 

3) A payload trim module which uses the same algo- 
rithm as the vehicle and determines a payload 
orientation which Is consistent with the vehicle 
trim flight conditions. 
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4) A nodule to calculate the payload state vector 
derivative which is used both by the trlnner and 
by the Integrator during the tine history* 

The combined trin is achieved by establishing the payload trin con- 
ditions consistent with the vehicle trim input flight conditions* The 
payload orientation is then calculated to meet the payload trim condi- 
tions* On completing this, the values of the cable forces are passed to 
and become part of the vehicle trim calculations* 

The vehicle trim is based on finding commands (6 linked commands) 
which will balance all forces acting on t'ae vehicle and result in no 
accelerations at the user input flight conditions* The payload trim 
uses the same algorithm, but attempts to find a position (linear and 
angular) where all forces are balanced* The accelerations will be con- 
sistent with the hull/payload rigid body motion determined by the vehi- 
cle trim flight conditions* 

The stability derivative matrices include both the vehicle and the 
payload states. A detailed discussion of the stability derivatives is 
in Part E of this section* The time history Integration loads the pay- 
load and vehicle state vectors into one vector and integrates them for- 
ward together. During the time history calculations, the cable spring 
and damping forces are the "elastic" constraints which control the pay- 
load position below the vehicle* 

B. VEHICLE-PAYLOAD IHPOT DATA FILES 

1* Data Piles: GMDIA, AEu*iXA, IFCDTA, FLMDIA, 

TBMDTA, and HISDTA 

These data files are used just as in program HLASIM. 

2* Data Files: PAYDTA 

PAYDTA contains all the basic payload data including: geometry, 

mass, aerodynamic and time history variables. (There are no inter- 
.ence effects on the payload). These inputs are listed in Appendix A 
a user will have to refer to the Technical Manual for the deriva- 
tion 01 ;he engineering constants. The form of the data very closely 
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paralltls Chat of Che vehicle and LPUs* The payload InpuCs are slallar 
Co Che hull, and Che Inpucs for Che fc .. cablee are similar In foraac to 
Che Inpucs for Che four LPUs* For mosc of Chese InpuCs, the descrip- 
tions given In Appendix A and In Che Technical Manual are sufficient* 
Only Inputs which have special considerations will be discussed here* 

Figure 12 illustrates Che four cable positions with respect to the 
payload and Che vehicle* Note, Chat "attach point 1" on the hull 
(RATHPl) oust correspond with "attach point 1" on the payload (ElPTCHl) 
and Chose attach points, by definition, are connected by "cable 1" (no 
crossed cables allowed)* 

The user may eliminate any of the four cables In a manner similar 
with the landing gears* To sero off any cable, the spring constant 
(CABLK) for that particular cable should be set to 0*0* The program 
will Chen calculate the force exerted by that particular cable to be 
zero, effectively removing It from Che model* The geometry Inputs for 
Chat cable must still be In the file but will not affect the calcula- 
tions* 

The restrictions on the payload cable geometry are: 

1) All the cable attach points on the payload may 

not all be on a common line segment Chat Inter- 

sects the payload center-of-gravity* This condi- 
tion cannot be trimmed. 

2) If the relative velocity in x-y plane or NUVP Is 
zero, all the active cables must not b*? attached 
Co the same point on the vehicle ot payload as 
this removes all yaw control and can not be 
trimmed. 

3) The user must Insure that none of the active 

cables (CABLK non-zero) will be slack at the trim 

position. The trimmer demands that all active 
cables have some tension In them* 

4) "Cable 1" must connect "Hull attach point 1" with 
"Payload attach point 1", and similarly for 
cables 2, 3 and 4 (l.e*. No crossed cables). 

These restrictions are Illustrated In Fig. 13. 
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Hull cable attach points (1-4) in hull 
reference axis 
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Figure 12, Payload Position Relative to Hull - Top View 







The cables (modeled as springs and dampers) cause Che trim 
calculation to be much more difficult Chan that for the vehicle. The 
user's choice of a cable geometry and spring stiffness and damping 

constants may sometimes cause the trimmer to fall to converge. These 
conditions are not predictable In that they arise from numerical 

problems Involved with the trim algorithm. There may be cases where the 
user may find it Impossible to model a certain geometry. Figures 14, IS 
and 16 includes Chose geometries which have been Implemented 

successfully during program development. 

The payload (l-mlnus-cosine) gusts are Input In Che same manner as 
those for Che vehicle (see namelist NPYGCOH). The payload gust strings 
are slightly different. The payload gust strings are turned on by Che 
flag PGSTFL similar Co the vehicle, but payload gust strings are read 

from two files. The user should read Che section on the Input files RG5 

and RG6 for a description of these files. The linear velocities are 
read from RG5, and angular velocities from RG6. These linear and angu- 
lar velocity vectors act at the payload aerodynamic reference center and 
are scaled by the variables PGVSCF (linear gusts), and POGSCF (angular 
gusts). 

The namelist NINDPST contains a set of variables which are used dif- 
ferently than the other Input variables. These variables, which are 
orientation perturbations away from the calculated payload trim condi- 
tion, are applied at the beginning of the t>oui history (time ** 0.0). 
These inputs allow the user to displace the p^iyToad from its trim condi- 
tion and study the resulting actions o2 the coupled vehicle/payload 
system during the time history. The payload perturbation variables are: 

DVPYLD - The displacement from the payload trim velo- 
city. 

DHRPYL - The dlsplacemenc. from the payload trim loca- 
tion. 

DPYELR - The displacement from the payload trim Euler 
rates. 

DPYEUL - The displacement from the payload trim Euler 
angles. 
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ORIGINAL PAGE 13 
OF POOR QUALITY 


$NRPTCH 

RPTCHl 

m 

40.0, -6.0, -6.0 

RPTCH2 

■ 

40.0, 6.0, —6.0 

RPTCH3 

m 

-40.0, -6.0, -6.0, 

RPTCH4 

$END 

$NRATHP 


—40.0, 6.0, —6.0, 

RATHPl 

- 

36.0, 0.0, 50.0, 

RATHP2 

- 

36.0, 0.0, 50.0, 

RATHP3 

m 

-36.0, 0.0, 50.0, 

RATHP4 - 
$END 

$NUSCLTH 

-36.0, 0.0, 50.0, 

USLTHl 

m 

22.0, 

USLTH2 

m 

22.0, 

USLTH3 

m 

18.0, 

USLTH4 » 
$END 

$NRPAYCG 

18.0, 

RPAYCG 

$END 


0.0, 0.0, 0.0, 

$NCABLK 

CABLKl 

m 

.9643E + 04 

CABLK2 

- 

.9643E + 04 

CABLK3 

$END 

$NCABLC 


.9643E + 04 

CABLCl 

- 

.3096E + 03 

CABLC2 

m 

.3096E + 03 

CABLC3 

- 

.3096E + 03 

CABLC4 

$END 

• 

.3096E + 03 


Figure 14. Cable configuration 1 (Inverted "V" and 
Front Lower then the Rear) 
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ORtGINAL PAGE B 
''F POOR QUALITY 


$NRPTCH 

RFTCHl 

m 

40.00, 

-6.00, 

-6.00 

RFTCH2 

m 

40.00, 

6.00, 

-6.00 

RFTCH3 

- 

-40.00, 

-6.00, 

-6.00 

RFTCH4 

$END 

m 

-40.00, 

6.00 

-6.00 

$NRATHP 

RATHFl 


36.00, 

-6.00, 

50.00 

RATHF2 

m 

36.00, 

6.00, 

50.00 

RATHF3 

- 

-36.00, 

-6.00, 

50.00 

RATHF4 

$END 

m 

-36.00, 

6.00, 

50.00 

$USLTH- 

USLTHl 


22.0000, 



USLTH2 

- 

22.0000, 



USLTH3 

- 

18.0000, 



USLTH4 

$END 

m 

18.0000, 



$NRPAYCG 




RPAYCG 

$END 


0.0, 0.0, 

0.0, 


$NCABLK 

CABLKl 

m 

.3100E+05, 



CABLK2 

» 

.3100E+05, 



CABLK3 

m 

.3100E-K)5, 



CABLK4 

$END 

m 

.3100E-K)5, 



$NCABLC 

CABLCl 

m 

.1243E+04, 



CABLC2 

m 

.1243E+04, 



CABLC3 

- 

.1243E+04, 



CABLC4 

$ENO 

m 

.1243E-K)4, 




> 


> 


Figure 15. Cable Configuration 2 (Four 'lables) 
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ORIGINAL PAGE 18 
OF POOR QUALITY 


i 


$NRPTCH 
RPTCHl - 
RPTCH2 - 
RPTCH3 - 
RPTCH4 - 
$END 

$NRATHP 
R4THP1 - 
RATHP2 - 
RATHP3 - 
RATHP4 - 
$END 

$inJSLTH 
USLTHl - 
USLTH2 - 
USLTH3 - 
USLTH4 « 
$END 

$NRPAYCG 
RPAYCG - 
$END 

I $NCABLK 

* CABLKl • 

CABLK2 - 
CABLK3 - 
CABLK4 - 
$END 

$NCABLC 
CABLCl - 
CABLC2 - 

! CABLC3 - 

] CABLC4 - 

I $END 

1 

i 

Figure 16. 


20 . 00 , 0 . 00 , - 8 . 00 , 

20 . 00 , 0 . 00 , - 8 . 00 , 

- 20 . 00 , 0 . 00 , - 8.00 

- 20 . 00 , 0 . 00 , - 8.00 


20.00, 0.00, 50.00, 

20.00, 0.00, 50 00, 

-20.00, 0.00, 50.00, 

-20.00, 0.00, 50.00, 


17.5000, 

17.5000, 

17.5000, 

17.5000, 


0 . 0 , 0 . 0 , 0 . 0 , 


.9643E404, 

0 . 0 , 

.9643E+04, 

0 . 0 , 


0.3096E+03, 

0 . 0 , 

0.3096E+03 

0 . 0 , 


Cable Configuration 3 (Bifilar Cables) 
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3. Data File: PTODTL 


Data file PYOUTL contains the two lists of code numbers for pa/load 
output variables* This data file is identical in format with data file 
OUTLST. Each list of positive numbers is terminated by zero or a nega- 
tive number. The first list contains the code numbers of the payload 
variables, and the second contains code numbers of the cable variables 
to be printed. The cable variables, like those of the LFUs, will be 
printed in groups of four, one value for each of the four cables under 
the variable name* 

4. Data Files: B6S MD) 166 

These two data files have the same format as RG1-RG4. Each line of 
the data file contains four numbers. The first number is the time and 
the following three numbers are the gust source vector for that time* 
Data file RG5 contains linear gust velocities at the payload aerodynamic 
reference center. RG6 Inputs contains angular gust velocities at the 
payload aerodynamic reference center. Payload gusts strings are dif- 
ferent than the vehicle gust strings in that they are applied directly 
to the payload, without spatial Interpolation. This requires angular 
velocities to be input, not calculated, hence the need for RG6. As 
before, interpolation is used to obtain the gust values at the appro- 
priate simulation time, and coordinate transformations convert the input 
data to the payload reference axes. 

The reference center for these gusts like the vehicle translates 
with the payload, but it does not rotate with the payload. Conse- 
quently, the gusts input in file RG5 and RG6 are in terms of the iner- 
tial reference axis, although their center i^ ‘ ranslatlng with the 
payload. The gust values being input at the vehicle sources are not 
Interpolated down to the payload nor are the payload gust values inter- 
polated up to the vehicle. 

5. Data File: INPUT (TAPES Interactive Questiona) 

The interactive questions for the payload program are the same as 
those for the main program (HLASIM). 
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C. PROGRAM HLAPAY OOTPDT 


The output listing for this program contains ell the input values 
for program HLASIM as well as the payload module. Appendix C at the 
back of this manual contains a sample output listing (Appendix B con- 
tains the corresponding input data files). 

Following the list of input values, the payload trim conditions are 
printed. The "Payload Trim Case Number" is a counter of the t xas being 
calculated (similar > HLASIM). The trim algorithm control is printed 
next, indicating the converged solution or the best guess achieved dur- 
ing tiie trim calculations. 

The payload trim algorithm is the same as that for HLASIM, and the 
meaning of "Iterations", "Restarts", and "Convergence Criteria" are the 
iame. The variable HBFYLC is the vector indicating the payload location 
relative to tha hull, and the three variables, PAY PHI, PAY THETA and 
PAY PSI are thf three Euler angles of the payload. These six numbers 
are the payload trim controls which are calculated during payload trim. 

During the payload trim calculation, there are two conditions 
encountered which are flagged as model errors. The position limit, 
HRPLFL, indicates that a guess for the payload position violated the 
condition that s^l active cables must have some tension for a valid 
trim. This particular variable is often flagged a few times during the 
trim calculation and does not indicate there is anything wrong. If the 
trim failed to converge and this variable is a large number, then the 
user geometry is probably such that one of the cables is carrying very 
little weight and trim guesses frequently causes it to be slack 
(restriction 2, Section 7-B2). This can oe remedied by altering the 
cable geometry slightly. The second error condition, indicated by 
counter PSNGMT, is a singular matrix. A large number for PSNGMT during 
a failed trim calculation often Indicates violation of or near violation 
of payload geometry restrictions 1 or 2. 
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Following this, the i>ayload variables are printed. The tlae Indi- 
cated here Is a negative number and functions In the same manner as the 
time does for the vehicle (it will be Che same time as the corresponding 
vehicle trim). The variables printed are those variables Che user has 
chosen through the input data file, PYOUTL. 

The vehicle trim data follows the payload data. This order indi- 
cates the order of the trim calculations. The user should consult 
Section 3 and A for a discussion of the vehicle trim. 

Program HLAPAY will generate a plotting file (PLOT) containing all 
the output data. The plot file has structure similar with the file gen- 
erated by program HLASIM. The post-processor PPLOTF (part of Che Heavy- 
Lift Airship Simulation package) will read this file and produce time 
history files compaClMe with NASA/Ames plotting software. This file Is 
described In detail In Section 3 -a^ and 9 of this manual. 

D. HLAPAY TRIM UNITATIOHS AND DIFFICULTIES 

The payload trim Is much more difficult than Che vehicle trim. This 
is caused primarily by Che springs that are modeled li the cableS' The 
user's able geometry has a large effect on the trim. If the trim will 
not converge foi a particular geometry, there are several things vdilch 
should be tried: 

1) Insure chat none of Che cable geometry restric- 
tions art violated (see Section 7-B), 

2) A reduction In the spring stiffness constants 
will often Improve the payload convergence char- 
acteristics. 

3. The user can alter his cable geometry slightly, 
as the program sometimes encounters numerical 
conditions which It can not handle and a slight 
change will avoid Chose problems. 

The discussion of data file PAYDTA (Section 7-B-2) contains a 11?* 
of cable geometries which have been Implemented successfully during pro- 
gram development 
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B. HLAPAT SIABILm D8BIVATIVBS 


The scabilicy derivatives for the prograa HLAPAT coatala the saoe 
basic matrices as are found in program HLASIM. Each of the four main 
matrices (A'-Matrix, B-Matrlz, B'*41atrlx and C Matrix) have 24 rows 
Instead of 12 (12 payload states are added to the 12 vehicle states). 
Bach of the auxiliary matrices contain 28 rows instead of 24. The four 
additional rows contain the cable force derivatives. 

The matrices generated are: 

A and A Auxiliary Matrices Generated by perturbing 
the 24 (12 vehicle and 12 payload states). 

B and B Auxiliary Matrices — Generated by perturbing 
the 27 vehicle control element effectors (same as 
program HLASIM). 

B' and B' Auxiliary Matrices -** Generated by perturb* 
ing each of the six linked vehicle controls (same as 
prograa HLASIM). 

C and C Auxiliary Matrices — Generated by perturbing 
the 48 (42 vehicle and 6 payload) gust velocities and 
velocity gradients). 

The same testing is dene for nonlinearities and those found are 
listed on the output, as in prograa HLASIM. Because of the symmetry of 
the vehicle and likely cable configurations, many of the cable derlva* 
tlves will be flagged as nonlinear when the payload position is per* 
turbed (l.e., pertubatlons right and left both Increase cable tension). 

There are two methods to obtain stability derivatives of tae payload 
alone 

1) Make the vehicle much larger (order of 100), so 
that the ratio of payload to vehicle mass Is 
small. The input variables Involved are: 

r 

HULVOL 

MASHUL 

IHULXX 

IHULYY 

1HUL7.Z 

t IHULXZ 
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The values for MASHUL and HULVOL aMSt be such 
that the rotors are able to support the differ** 
ence between g(HASHUL) and the buoyancy. 

The resulting aatrlx rows and coluons pertaining 
to the payload only will be the decoupled payload 
stability derivatives. Slnilarly, the eigen*- 
values and eigenvectors pertaining to the payload 
are chose tor the decoupled payload systea only. 

2) Subroutine PLIMAR calculates the payload only 
stability derivatives. This sodule is in the 
prograu library, but it is not called in the 
overlay structure at present. The subroutine 
TLINAR, called froa the aain program HLAPAT, cal- 
culates the coabined vehicle/payload stability 
derivatives. If the user wishes to have only the 
payload stability derivatives calcula**ed, he 
should change this call to subroutine PLINAR. 

There are nc other changes necessary to the code, 
and the program will run as usual. The trim will 
be calculated for the payload and the vehicle, 
but stability derivatives will be calculated for 
the payload alone as if attached to an infinite 
mass vehicle. 

The above methods proviia information on the payload alone charac- 
teristics, without the complex coupling between the vehicle and payload 
systems. If a m-.w payload model is developed, these payload only sta- 
bility derivatives may prove useful in debugging and verifying the pay- 
load model. This would be much cheaper since the program will calculate 
only the elements of the payload only derivative matrices. The user 
should refer to the Programmer's Manual for detailed description of the 
structure of the stability derivative modules if he is contemplating 
making these changes. 
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SBCTIOH 8 


PIOGBAM HLAICA; SDflJIATION OF 1HK MD08BD OaFOiBIED 
HKAVT LIFT AHSHIP 


A. IHIKODDCTION 

This prcgraik simulates an unpowered heavy lift airship which Is 
translatlonally moored to a mast. Since the mooring point Is modeled as 
a transL'tlonally constrained perfect glmbal, the airship has only three 
degrees of freedom (angular glmbal motion). No linear motion at the 
mast attach point Is allowed. The vehicle Is unpowered and consequently 
has only minimal rotor and propeller aerodynamics. The trim Is con- 
siderably different ' than that for HLASIM, as It Is based on the mooring 
point constraint, the active landing gear forces, and the wind forces. 
There Is no flight control system and no control commands. The only 
data which disturbs the trimmed vehicle are the gusts and the mooring 
state displacements. 

B. MOORED VEHICLE INPUT DATA FILES 

This program uses the same basic data and consequently will model 
th** same vehicle as program HLASIM. Many of the Input data for program 
HLASIM are not used in the mooring program because of the unpowered con- 
dition. Even so, the Input files have been structured In such a manner 
that the same files that were used by program HLASIM may be used here. 
This program will Ignore those Inputs which are not relevant to It. 
This has been done to simplify the user's problem of Insuring the data 
used In HLASIM and HLAMOR are the same. 

1. Data Flies: GHDTA, ABOOTA, IFCDTA, 

TSMDTA, HISOTA 

These files may be used just as they are from program HLASIM. Note, 
file PLMOTA Is not used at all In this program. 


TR-1151-2-III 


67 



la file TRNDTA, the only necessary data is contained In namelists 
NATMOS and NSTABDV. The data In namelist NINSTAT Is not used. For this 
reason the flag "SUBROUTINE INATMOS" has been Inserted In this file Just 
proceeding the namelist "NATMOS". This allows this program to read past 
the unwanted data. It is highly recommended that the user, if he Is 
modeling the same vehicle In program HLASIM and HLAMOR, use the complete 
data file and not remove the unnecessary data. Otherwise, the user may 
Inadvertently alter data In the file, and Invalidate the comparison 
between different runs. 

Similarly In file HISDTA, the flight control system data Is not 
necessary to the program. The program reads past all the data until It 
encounters the flag "SUBROUTINE INGUST" indicating the beginning of the 
gust inputs. 

2. Data File: MORDTA 

This data file contains the variables which are exclusive to HLAMOR. 
The variables are listed In the tables In Appendln A. They are all In 
namelist format. 

The variable PSIO serves two functions: 

1) The vehicle Inertial Euler yaw angle is set to 
PSIO If there Is no wind. 

2) The vehicle Euler yaw angle relative to the wind 
for the trim Initial guess is set to PSIO when 
there Is any x*-y plane wind velocity. 

The moored vehicle trim may have three converged yaw positions: nose 
Into the wind or yawed positive or negative relative to the wind. The 
trimmer can be made to converge to any of the three positions by using 
PSIO to set the initial trim guess close to the trim yaw position 
wanted. PSIO Is the wind relative vehicle trim Initial guess Euler yaw 
angle if the x-y plane wind Is nonzero. PSIO « 0.0 causes nose Into the 
wind trim initial guess. The other two positions will depend on the 
user vehicle characteristics. 
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The trim solution for a no-vind condition is achieved by creating an 
artificial yawing moment* It is scaled to the difference between the 
desired yaw angle (PSIO) and the actual Euler yaw angle, and will only 
exist during trim (MFORCE, CLMTRM). This places the following restric- 
tion on the possible landing gear geometry. 

If there is only one landing gear, or if all active landing gears 
are coplanar and aligned parallel to the mast, then the composite vehi- 
cle center of gravity must lie in the (landing gear) plane* Otherwise 
the artificial moment generator will cause the trimmer to converge on a 
non-existent solution* This geometry restriction will only apply in the 
no-wind condition* 

The three variables DELTAL, ISLTEL, and l^LTRD are the tail deflec- 
tion angles for the ailerons, the elevator, and the rudder* These 
values are fixed positions which arc held throughout the time history* 

The variable OHLEUL is a perturbation vector away from the trim 
value of the Euler angles (similar to the variables in namelist NI>*I>i-ST 
of data file PAYDTA)* Its purpose is to allow the user to move the 
vehicle away from the trim condition at the beginning of the time his- 
tory and study the subsequent motion of the vehicle* Aside from input- 
Ing a gust, this is the only means of causing the vehicle to move away 
from the trim (during the time history)* 

3* Data Pile: ODTLST 

Program HLAMOR, like HLASIM, uses OUTLET to specify the output vari- 
ables to be printed* If the user requests variables which have not been 
calculated, the program will print either zero or an undefined ("I")* 
Because HLAMOR is unpowered and does not use the sophisticated aero- 
dynamic models for the rotors and propellers, mar.'' of the variables 
which can be requested are not being calculated* 

4* Data Pllefi: RGl, RG2, RG3, RG4 

These four gust string input files are used just as in program 
P'JiSIM* 
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PROGRAM HLAMOR OOTTOT 


The HLAMOR output listlog Is the sene as the HLASIM listing, except 
the input variables not used in HLAMOR do not appear. Output variables 
which are not calculated in this program should not be requested through 
file OUTLST. If they are requested they will appear as zero or unde- 
fined valuei>. 

The plotting file (PLOT) is the same as the HLASIM plotting file. 

D. HLAMOR TRIM OOHSIDEIATIOIIS 


The trim calculation for HLAMOR is a three degree-of -freedom prob- 
lem. The vehicle linear position is determined by the user input mast 
location and mooring point. Consequently, the vehicle can only move in 
the three angular directions. The trim is much easier and usually 
requires fewer Iterations than program HLASIM or HLAPAY. There are a 
few restrictions however. 

1) The geometry of the mooring point, mast location, 
and landing gears must be such chat all active 
landing gears will be touching the ground. The 
landing gear :rames, vehicle tail or belly, must 
not be touching the ground. Deactivating all 
landing gears will cause a buoyant (airborne) 
trim to be calculated. In this case Che buoyancy 
and wind must be enough to hold the vehicle belly 
off the ground. 

2) The only physical means for establishing a yaw 
angle is the force of Che wind acting on the 
vehicle. Consequently, the user must input a 
value (PSIO) indicating the yaw angle In the 
event that the wind vector Is equal to zero. 

3) If only one landing gear is active it must be in 
the same lateral position as Che composite center 
of gravity of the entire vehicle (hull and 
LPU's). 

E. HLAMOR STABILITY OBRIVATIVES 


Stability derl.atlves for the moored vehicle are reduced to the 
A-Matrlx and the C-Matrlx and their corresponding auxiliary matrices. 
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• A-^trlx — This 6 by 6 matrix is generated by 
perturbing the 3 Euler angles and the 3 Euler 
angle rates of the vehicle* 

• C-Matrlx — This 42 by 6 matrix Is generated by 
perturbing all of the gusts acting at the various 
component positions (hull, tail, LPUl-4); Identi- 
cal to HLASIM. 

The auxiliary matrices each contain 28 rows, consisting of the 24 
constraint forces and moments between the LPUs and the vehicle, and the 
four compression forces acting through the landing gears. 
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SBCTION 9 


FLOTTIHG FILE FOST-FSOCSSSOt FIDGIAM FPLOIP 


Prograa PPLOTF (process the plotting tile) 'trill reforaat the plot- 
ting file (PLOT) produced by any of the slmulatica programs (HLA31M, 
HLAMOR, or HLAPAY). The Input for program PPLOTF Is the binary format 
file PLOT* The output Is a binary file (THPLOT) of the tl»ie history 
data In a format con^atlble with the NASA/Ames flight data plotting 
software. The Heavy-Llft Airship Simulation data will be plotted with 
the NASA/Ames flight data plotting software. 

PPLOTF program trrltes to the standard output flic, OUTPUT, a title 
which Includes the simulation name (HLASXM, HLAMOR or HLAPAY), the date 
of the simulation and the simulation times (Trim and Time history). 
This Informs the user which simulation run plotting file Is being pro- 
cessed* If this program finds something wrong with the Input file 
(PLOT) It will print a message and stop. 

File PLOT contains all the trim data, variable names and time hls- 
tor^r data. Program PPLOTF Ignores the trim data. The variable names 
contained in file PLOT are exactly as they appear on the simulation pro- 
grams' output listing. They are changed slightly when written to file 
THPLOT. 

1) Each LPU and cable variable name has four values 
(LPUl-4 and Cable 1-4). This program writes each 
name four times Inserting numbers 1 to 4 In the 
blank sixth position. 

For example RFIV :X becomes 

RFIV 1:X, RFIV 2;X, RFIV 3:X, and RFIV 4:X. 

2) All leading and embedded blanks are then removed. 

RFIV becomes finally: 

RFIV1;X, RFIV2:X, RFIV3:X, and RFIV4:X 
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3) The independent variable. Time, is given the name 
"ETIME" and ^ts value la written with the other 
variables in each data frame. 

These alterations are not seen by the user until particular varia- 
bles are plotted with the NASA/ Ames flight data plotting software. 

This program is executed by loading the program with the data file 
PLOT. The file THPLOT should be saved after execution is complete. 

Secti%;a 3-B contains a detailed description of the file PLOT. The 
flJi! THPLOT is outlined in the diagram below. 


Record 1 Time step, Number of variables. Date (Julian) 

and Variable names 


Record 2 Data time frame 1 

Record 3 Data time frame 2 


(Continued until end of time history) 


Note: The number of variables includes the time, and the 

variable names Include "ETIME" in the appropriate 
place. The time is written as a data value in each 
time frame. 
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ORIGINAL PAGS 13 
OF POOR QUALITY 


SBCTIOH 10 

LINBARlZmOH NmiX ROM GOST SOOICBS 10 
VEHICLE OQMPONEHTS 


Program GSRCSB is provided to aid gust model analyses* It generates 
a gust source interpolation matrix (GSI) which shows the relationship 
between the gust velocities at the four gust sources and the gust velo- 
cities and gradients at the vehicle components. Using the vehicle and 
gust source geometry and che trim states, the program perturbs each gust 
source velocity and loads the resulting vehicle component gust veloci- 
ties and gradient derivatives into columns of the GSI matrix. The fol- 
lowing equation defines the GSI matrix. 
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( 10 - 1 ) 


where DUGBAR Is the same vector of gust velocities, accelerations, and 
gradients found in the main program C-matrlx stability derivative calcu- 
lations. 


This matrix will change only when the vehicle geometry, gust source 
geometry or hull Euler triu angles are changed. Consequently, program 
GSRCSB iihould be executed only if some of those variables have been 
altered. 

The program uses the main simulation Input routines and can read 
the main program data files. The files required are: GMDTA, ARODTA, 
TRMDTA, HISDTA and ERMSSG. The output is written to file OUTPUT. 
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Program GSRCSB does not calculate the hull or tall gtist acceleration 
rows (U DOT, V DOT, •••, R DOT) of 6SI. For example, thu row corres- 
ponding to the hull gust acceleration (u^j term Is: 

GSI(1,1) * D/DT, GSI(1,2) * D/DT, •••, GSI(1,12) D/DT 


'*'nis nomenclature symbolizes the following mathematical expression for 
u^ denotes row 1, column 1 of the GSI matrix): 


• _ du®^ . dv®^ . dw®^ . du®2 

uh “ GSIi,2 GSIi,3 GSIi,4 


dv®^ dw®^ 

+ — + GSI1.11 _+GSIi,i2-^ 


(10-2) 


and similarly for the remaining hull and tall gust accelerations. 


The main prog .m stability derivative C-matrix Is generated accord- 
ing to the following equation (see stability derivative output Sec. 5). 
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(10-3) 


The GSI matrix Is used with the main program linearization C-matrlx 
In the following manner. Equations 10-1 and 10-3 are conceptually com- 
bined to form £q. 10-4, which shows the effect of gust source velocities 
on the vehicle. 
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(10-4) 
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The method to mealpulate the d/d^ terms of the gust scceleratloa 
rows (GSI matrix) in the calculation of Eq. 10-4 will depend on the 
s''ftware utilities available to the user. 
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ovEuufs, EZBorrioH or the ncraEis w ibb iasa/aibs 
CDC 7600/S00FB SYSTEM 


The program has been Installed on the NASA/Ames CDC 7600/cCOPE sys- 
tem in the following manner: 

1) Simulation program source code Is stored on disk 
using the UPDATE facility. 

2) All data files are stored as separate disk files 
with the names given In Section 2. 

3) The three main program routines (HLAMOR, HLASIM, 
and HLAPAY) are stored as separate disc files as 
well as In the UPDATE facility (see Item 1). 

4) The segmented load directive files are stored as 
separate disc files (SLDIR, for HLASIM; SLDIRM, 
for HLAMOR, and SLDIRP, for HLAPAY). 

5) The binary object code for all program subrou- 
tines (excluding main program routines) Is main- 
tained In a Library on a a disk file called 
MAINLB. It can be accessed and edited with the 
LIBRARY facility. 

6) The post processing programs (<*PLOTF and 6SRCSB) 
rre in separate disc files. The files are named 
PPLOTF and GSRCSB, respectively. 

A. PROGRAM EZECUTIOIf OCWTROL SEQOENCBS 


Figure 17 (page 79} Is a sample job control sequence for executing 
the main simulation ,/rogram. A line by line discussion of this sequence 
follows: 

1-2 Job control and accounting information. 

3 The IMSL library le attached. 

4-5 The disc containing the program files Is Identified and mounted. 

6 The plotting file will be written here (If It is generated). 
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7^23 Data files are attached. J.1 data files beiog attached are from 
the basic set (see Section 2). If a special geoaetry was to be 
inpleittcnted a copy of GMDTA Mc^xld be aad*^ containlag the desired 
change. This file» naned GHliXTA, would be attached by: 
ATTACU(GMDTA,GM1DTA» ID»userID) . 

Similar changes could be introduced in the other data files. To 
ease the "bookkeeping” problems, it is recoisfiended "hat the 
notation Just introduced be continued (naise other data files 
AROIDTA, etc.). In this case the output fil ' could be called 
RUNl. 

25 The filename in this line should be HLASIi., HLAMOR, or HLAPAY. 

26 The main program routl'ie attached in line 25 is conpiled. 

27 LMAP is a file for the load map. It will be saved only if the 

program aborts. 

28 MAINLB and IM£L, Libraries for usa by the loader, aro ciclared. 

29 Files ricntaining the segmented load directives attached. 

The file attached xmx^t correspond with the i-togi . o be run 
(SLDIR, for KIASIM; SUIFM, tor HLAMOR; end SLDIRP for HLAPAY). 

30-34 )*uad and execure sequence. 

35 Normal termination of Job. 

36-37 In the event of abnormal Job termlnatioi: dxx^v tl»c load map onto 
output file. 

39-49 Resporses to interactive questiouc (see Sec. 2-1;. 

Between lines 34 and 35 

If pl^ ting files are being generated the following line should 
be inserted at this point. 

CATALOG ( PLOT , PLOT 1 , ID»use r ID ) 

This saves the plotting file as PLOTl (continuing the notation 
from above). 

Alternatively, the complete command string could be inserted 
here to ''xeente PPLOTF (see below). 

The plotting file post processor, PPLOTF, can be erecuted with the 
Job control sequence shown in Fig. 18. The filename in line 6 is the 
plotting: file written by the main simulation progrnu: Following the 

c*. ample given above, the filename would be PeOTl. THPLOT is the output 
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file and is cataloged in line 10. Following our notation, the fllenane 
would be THPLOTl. This file is read into NASA/Anes flight data plotting 
software. 

The gust source stability derivatives prograa can be executed using 
the control sequence given in Fig. 18 with a few changes. 

1. Make filename in line 2b to be GSRCSB. 

2. Delete line 31. 

Many of the other lines are unnecessary. If the user wishes to delete 
those lines. Section 10 of this manual should be consulted. 

B. SUaiEllTBD LOAD DIBBCTIVKS 

The segmented load directives for each of the three programs are 
listed in Figs. 19, 20, and 21. They are listed here without explana- 
tion; they are discussed in the Programmer's Manual . 

As noted earlier, the appropriate segmented load directive file must 
be loaded with each of the three programs. 
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AIRSHIP. T30O ,PN , 70 1 , YL 1 
ACCOUNT , user ID» T 1498L2 
ATTACH (IMSL . IMSLLIB , ID^^AMESLIB ) 
SETNAME(HHLASqR) 

MOUNT (VSN-D0413A) 

REQUEST (PLOT, *PF) 

ATTACH (PYOUTL, PYOOTL, ID»userID) 
ATTACH (PAYOTA, PAYDTA, ID«user ID) 
ATTACH (GMDTA,GMDTA, ID«userID) 
ATTACH (ARODTA,ARODTA, ID-userlD) 
ATTACH (HISDTA, HISDTA, ID-user ID ) 
ATTACH (PLMDTA, PLMDTA, ID-user ID) 
ATTACH (TRHDTA, TRMDTA , ID-user ID ) 
ATTACH ( IPCDTA, IPCDTA, ID-user ID ) 
ATTACH (MORDTA,MORDTA, ID-user ID ) 
ATTACH (ERMSSG, ERNSSG , ID-user ID ) 
ATTACH (OUTLST , OUTLST , ID-user ID ) 
ATTACH (RG 1 , RG 1 , ID-user ID ) 

ATTACH (RG2 , RG2 , ID-user ID ) 

ATTACH (RG3 , RG3 , ID-user ID ) 

ATTACH (RG4,RG4, lO-userlO) 

ATTACH (RG5 ,RG5 , TD-user ID ) 

ATTACH (RG6,RG6, ID-userlD) 
ATTACH(MAINLB ,MAINLB , ID-user ID ) 
ATTACH (PG , filename , ID-user ID ) 

FTN (I-PG ,L-0 . PL-5aM0) 

REQUEST (LHAP,*PF) 

LIBRARY (MAINLB , IMSL , *) 

ATTACH (LOADIR, filename , ID-user ID ) 
MAP (ON) 

SEGLOAD(I-LOADIR) 
LDSET,MAP»X/Ll«AP, PRESET-INDEF. 
LOADaCO) 

EXECUTE 

EXIT 

REWIND (LMAP) 

COPY (LMAP, OUTPUT) 

ZZEOR 

T 

01 


T 

T 

SIX LINES OF COMMENTS 
TO APPEAR AS RON DESCRIPTION 
ON OUTPUT LISTING 


Figure 17. Job Control Sequence to Load and Execute Program 
TR-1151-2 tit 80 


• - ( 



OHJGINAL PAGE fS 
OF POOR QUALITY 


1 AIRSHIP, T50,PN,YD1,YL1 

2 ACCOUNT, PVSHM.T1498L2 

3 SETNAME(HHLABQR) 

4 MOUNT(VSN-D0419A) 

5 ATTACH(PPLOTF,PPLOTF,ID«FVSHM) 

6 ATTACH(PLOT, filename, ID-FVSHM) 

7 REQUEST (THPLOT,*PF) 

8 FTN(I»PPL0TF,L»0) 

LGO 

10 CATALOG(THPLOT, filename, ID«FVSHM) 

11 EXIT 


Figure 18. Job Control Sequence to Execute 
Program PPLOTF 
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ORIGINAL PAGE 15 
OF POOR QUALITY 
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Figure 19. Segmented Load Directive File (SLDIR) for Program HLASIM 
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Figure 20, Segmented Load Directive File (SLDIRM) for Program HLAMOR 
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Figure 21, Segmented Load Directive File (SLDIRP) fcr Program HLAPAY 



